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Abstract 

In this article we review recent theoretical and experimental devel- 
opments on multilongitudinal-mode emission in ring cavity lasers, pay- 
ing special attention to class B lasers. We consider both homogeneously 
and inhomogeneously broadened amplifying media as well as the limits of 
small and large cavity losses {i.e., we treat cases within and outside the 
uniform field limit approximation). In particular we discuss up to what ex- 
tent the experimental observations of self-mode locking in erbium-doped 
fiber lasers carried out in recent years are a manifestation of the Risken- 
Nummedal-Graham-Haken instability. 

1 INTRODUCTION: 

WHY SHOULD WE CARE ABOUT LASER 
INSTABILITIES? 

The very first laser - Maiman's ruby in 1961 - produced extremely unstable 
emission, as evidenced by a figure in the original publication (Ij. Soon there- 
after, researchers began to learn the skills of how to avoid instabilities in lasers 
designed for applications. Nonetheless laser instabilities have been around ever 
since both as a nuisance lurking to haunt technically-minded people, and as an 



interesting nontrivial physical phenomenon for the more fundamental-principles- 
minded set. After all, lasers are nonlinear dynamical feedback systems, and 
instabilities are inherent in such systems. It is certainly warranted to gain an 
understanding of laser instabilities - if only for the minimalist purpose that 
they be avoided successfully in spite of the ever increasing demands on lasers in 
terms of power, speed, tunability, etc. 

In lasers, typically one or several modes of the light field are subject to res- 
onator boundary conditions while at the same time being dynamically coupled 
to the amplifying medium. Typically the coupling is highly nonlinear, and de- 
pending on the laser type and particular circumstances, the laser may behave 
in many different ways, running the gamut from smooth and stable single- mode 
to irregular and unpredictable chaotic operation. 

To systematically address the plethora of possibilities, one can make a first 
distinction between cases in which either several transverse modes are involved, 
or just a single one (TEMoo). Next-neighbor longitudinal modes always have 
very nearly the same frequency difference, and only a single beat frequency, 
along with its overtones, can occur. In contrast, for transverse modes there can 
be more beat frequencies, giving rise to an enormous wealth of possibilities for 
frequency locking and pulling phenomena. Maiman's ruby laser, for example, 
displayed an instability that involved several transverse modes. The complexity 
of the situation is not helped by the fact that a mathematical description requires 
either an infinite set of ordinary differential equations, or a set of fully space- 
dependent partial differential equations. Nonetheless, quite some research was 
devoted to multi-transverse mode dynamics (see |21E1E] for references, and the 
rest of articles appearing in the present volume). More attention, however, was 
given to single mode laser instabilities over the last three decades. Overviews 
and references can be found, e.g., in pi 1^ 151 1^ [71151 1^ [TUl [TTl [HI [T^ . 

Here we shall consider single-transverse mode problems. In this context, 
there is an important distinction whether several longitudinal cavity modes are 
involved in the lasing process, or just a single one, although a connection between 
singlemode and multimode instabilities can be established It is also of 

fundamental importance whether the gain medium exhibits a homogeneously or 
inhomogeneously broadened lasing transition (see, e.g., [21 El)- 

It was noticed by Haken in 1975 jJS] that the MaxwcU-Bloch model for a 
homogeneously broadened single mode laser is isomorphic to the Lorenz model 
of chaos, originally developed for climatic instability ^H]- This directly implies 
that these lasers must exhibit a second threshold, i.e. a characteristic value of 
the pump power above which a Lorenz-type instability sets on. One can show 
that this second threshold is at least nine times higher than the first (ordinary) 
threshold above which there is coherent oscillation (see, e.g., pil^ IIUI 1151 
This "factor-of-nine" is of great importance to the present paper. Later Weiss 
found ^7j that certain far-infrared lasers have suitable damping rates and can 
be pumped hard enough to test this prediction. Indeed they displayed chaotic 
behavior remarkably similar in many respects to that exhibited by the Lorenz 
model Unfortunately the level structure of the gain atoms was much 

more complex than a two-level system. The relevance of this difference for 
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the laser dynamics was the focus of some controversy during the eighties and 
the beginning of the nineties (see [201 1^ for a discussion at depth and for 
references). However, in the absence of a formal proof of equivalence between 
the realistic laser model and the Lorenz model, one can only cautiously conclude 
that some lasers can display Lorenz chaos even when their structure does not 
make it obvious that the Lorenz equations are the appropriate model. 

In contrast, multilongitudinal mode laser emission has received only marginal 
attention from the viewpoint of laser instabilities. This is all the more remark- 
able considering that most lasers emit in several longitudinal modes, if only due 
to spectral hole burning in inhomogencously broadened active media (see, e.g., 

1221113 El) 

It is a common belief that with a homogenously broadened gain line and 
single-transverse mode condition, the only way to have more than one longi- 
tudinal mode oscillating is spatial inhomogeneity in the medium, such as the 
spatial hole burning occurring in linear resonators. The spatial inhomogeneity 
provides the required amount of independence of the inversion available to one 
mode from that available to another mode so that the usual winner-takes-all 
coupling is suspended in favor of the mild mutual coupling due the finite overall 
energy balance^. 

This myth was shattered when in 1968 two publications, independent of 
each other, discussed the situation in detail |2HI IHOI (see also 123 )• It was 
shown that even in a unidirectional resonator filled with a perfectly homogeneous 
medium, Rabi splitting of the lasing transition induced by the lasing mode can 
provide gain for other longitudinal modes. This means that even perfect gain 
homogeneity does not safeguard against instability. While a string of four names 
makes a somewhat awkward moniker, fairness dictates to call this mechanism 
the Risken-Nummedal- Graham- Haken Instability, or RNGHI for short. Ikeda et 
al. introduced the term Resonant Rabi Instability for designating this instability 

In 1201 12] a unidirectional ring cavity with small cavity losses filled 
with a homogeneously broadened two-level active medium was assumed. The 
prediction was that for RNGHI to occur the pump power must exceed a certain 
instability threshold. This instability threshold is referred to as the second 
threshold, in distinction to the familiar first threshold which defines the onset 
of coherent laser oscillation. The value of the second threshold came out to be 

good deal of theoretical and experimental work has been done in order to correctly 
understand multimode emission in Fabry— Perot cavity lasers, starting with the well known 
Tang— Statz-de Mars model 11^51 . We refer the interested reader to ll2ll^till^7ll^SI and references 
therein. That subject, while very interesting in itself, is beyond the scope of the present 
discussion. 

^Interestingly, in 1976 Graham 1331 showed that the multimode laser model, the RNGH 
model, is isomorphous to the Lorenz model |16|. thus extending the analogy discovered by 
Haken the previous year 1151 . The difference between the singlemode and the multimode case 
is that in the latter the parameter a is not fixed as it depends on the velocity of the travelling- 
wave solution. This would be extended further to cover detuning in 1990 when Ning and Haken 
[341 showed that the detuned multimode laser equations are isomorphic to the complex Lorenz 
model for the baroclinic instability 35 . This isomorphism is a powerful tool that has not 
been investigated enough, specially for what concerns the multimode instability. 
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at least nine times the lasing threshold (the same "factor-of-nine" as in single- 
mode instabilities in the Lorenz-Haken model). Another condition was that 
the laser cavity length must exceed a certain minimum. Indeed, it has to be 
unrealistically long for conventional bulk lasers |21 E] . 

Due to these predictions and without experimental demonstration of the 
opposite, interest in RNGHI abated over time: While during the Seventies and 
Eighties much theoretical effort was devoted to the understanding of the RNGHI 
as well as to the dynamics of the laser above the second threshold [201 El UHl 
EHl EHl ^1 im J eventually RNGHI was more or less dismissed as a merely 
academical prediction rather than an actual mechanism for multimode emission. 

In 1984 researchers announced that they had observed what could be inter- 
preted as RNGHI in a dye laser ^B] . The dye laser seemed to be a good 
candidate for the observation of the RNGHI as (i) it is a homogeneously broad- 
ened laser, and (ii) the "critical" cavity length (see later for a precise definition 
of this term) is smaller than typical resonator lengths 0. However, it was pecu- 
liar that in the spectrum there where only two peaks instead of three. Finally, 
after some debate 0^1 ; Fu and Haken showed that a suitable model for the dye 
(consisting of bands rather than levels) could explain the observations without 
invoking the RNGHI gTlEHlEni- From then on the observation in [21113 was 
no longer considered as a manifestation of the RNGHI. 

A few years later, a phenomenon very similar to RNGHI was discussed in 
the context of optical bistability performed in the microwave regime |5(JI 151) . 
The experimental system consisted of a waveguide Fabry-Perot cavity with a 
length of 182 m filled with hydrocyanic acid vapor at low (« 1 mTorr) pressure. 
Driven at a frequency of 86 GHz near a HC^^N transition, this system displayed 
bistability and, on the upper branch of the bistable loop, self-oscillation at a 
frequency of the order of the cavity free spectral range and of the Rabi frequency, 
but could deviate by 50% or so. Nevertheless it is to be emphasized that this 
fairly unique experiment actually displayed not the RNGHI, but the multimode 
instability of optical bistability predicted in 1978 by Bonifacio and Lugiato 
Thus, for many years and after a false start, the closest thing to the RNGH laser 
instability that was actually observed was Rabi splitting in microwave optical 
bistability - not very close, certainly. 

The discussion on the RNGHI during the late Eighties and the Nineties 
faded, although theoretical studies continued pilHHllK^IHHllK()llF7llKHllK??llHni 
inilE2ESlEllEHllMlinZlEH|,and signatures of the RNGHI were predicted also 
for laser related systems such as a laser with an intracavity parametric amplifier 
1^. and a laser with injected signal [7D]. Nonetheless, in this period it seemed 
like the matter of longitudinal mode instabilities in lasers, and the RNGHI in 
particular, would soon be more or less forgotten. 

The situation changed with a suggestion by Lugiato and coworkers |71l I72| . 
They observed that a ring-cavity erbium-doped fibre laser (EDFL) sponta- 
neously mode-locked and emitted a train of pulses with a repetition time equal 
to the cavity roundtrip time. This is one signature of the RNGHI, but certainly 
not in itself a sufficient criterion. However, the unusually long cavity length 
of fibre lasers automatically fulfils the most difficult prerequisite for RNGHI of 
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sufficient resonator length. On this basis they put forward the hypothesis that 
self-pulsing in these lasers could be a manifestation of the RNGHI (see also [75] 
for compatible experimental results). 

Unfortunately they could not confirm the existence of a second laser thresh- 
old, i.e. a threshold-like onset of the instability. Moreover, there was a grave 
quantitative difficulty: Instabilities were observed at pump powers immediately 
above the lasing threshold [72], in stark contrast to the prediction that the 
second threshold be at least nine times higher than the first. 

It was subsequently shown that this requirement is not necessarily applicable 
to real-world lasers. In [TJ] the usual two-level atom model for Lorenz-type 
instability was extended to incorporate a third level as a necessity for applying 
a pump source. It turned out that the instability threshold was lowered such 
as to come close to the lasing threshold. This finding was then applied to 
EDFLs [231 • Erbium ions, too, behave more like three-level systems [JH]; rather 
than two-level systems as the original theory assumed. Now the ratio of the 
instability to the laser threshold was predicted to be close to unity. This was 
subsequently corroborated (and expanded to include four- level cases) in |77) . 

Put simply, for EDFLs the expected "factor-of-nine" , became a "factor- 
of-(l -I- e)". Thereby, the apparent contradiction with the hypothesis in [72] is 
gone. This insight bestowed fresh vigor on the debate on RNGHI and lead to 
dedicated experimental investigation as will be described below. 

In this paper we shall review the fundamentals of the RNGH theory as 
well as the research we have carried out along the recent years (including some 
inedit results), i.e., we shall not review in detail the available literature. After 
this historical introduction, the paper follows with three more sections. Sect. [2] 
is devoted to the RNGHI in homogeneously broadened media. In this section 
we first explain in detail modelling issues. Sect. 12.11 and then we treat the 
RNGHI, first in the uniform field limit (Sect. I2.2|l and then outside the limit 
of validity of this standard approximation (Sect. I2.3|) . In Sect. [3] we consider 
multilongitudinal mode emission in inhomogeneously broadened media, again 
within (Sect. I3.2|l and outside (Sect. the uniform field limit. Then in 

Sect. 0]we discuss the experimental aspects of the RNGHI concentrating our 
discussion on erbium-doped fibre lasers. Finally, in Sect. [S]we provide a general 
discussion of the issue. 

2 MULTILONGITUDINAL MODE EMISSION 
IN HOMOGENEOUSLY BROADENED RING 
LASERS: THE RNGH INSTABILITY 

2.1 Modelling 

In this section we introduce the standard two-level laser theory (Sect. 12.1.1(1 . 
We then discuss how the model can be applied to three- and four-level lasers 
(Sect. 12.1. 2(1 . and how it can be generahzed to treat the dependence of the 
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laser field and of the pump field on the radial coordinate (Sect. I2.1.3|l . Finally, 
we determine the stationary singlemode solutions of the laser equations (Sect. 
12.1.41 and we rigorously derive the laser equations in the uniform field limit 
(Sect. [TT3|l . 



2.1.1 Two level atoms 

We consider a collection of N two-level homogeneously broadened atoms per 
unit volume, with transition frequency Wa interacting with a linearly polarized, 
plane wave, unidirectional electric field of carrier frequency Wc, chosen as the 
longitudinal mode frequency closest to to a- The Maxwell-Bloch equations de- 
scribing the system are [TTl [T21 gSl EH] 

d.F + v:^'d,F = ^P-^F, (1) 

dtP = 7i [FD - (1 + i5)P] , (2) 
dtD = [1 - D - Re {F* P)] . (3) 

Here F, P and D are properly scaled variables representing the electric field, 
medium polarization, and population difference 

F=-^, P--2.^,/^, D^^. (4) 
VTlTiT ao V 7|| "0 

In these expressions E is the Rabi frequency associated with the laser field, pi2 
is the slowly varying envelope of the coherence between the lasing levels, and d 
is the population difference per atom between the upper (2) and lower (1) lasing 
levels; 7j_ and 7|| are the decay rates of the medium polarization and population 
difference, respectively, S = (wa — Wc)/7j_ is the detuning between the atoms and 
the cavity, Vm = c/nm is the light velocity in the host medium with refractive 
index n^, am is an intensity loss coefficient per unit length, which describes 
(non resonant) distributed losses inside the active medium |78| . and c?o is the 
equilibrium value towards which d relaxes in the absence of an electric field, do 
is positive because the medium is amplifying. The unsaturated intensity gain 
coefficient per unit length a is proportional to the total population difference 
Afdo, and we can write it as 

a = — TV do , (5) 

where /z is the dipole moment between the two levels of the laser transition. 

Equations must be supplied with the appropriate boundary condition 

for the electric field. We assume that the active medium fills a region of length 
Lm inside a ring cavity of length Lc- Denoting by z = and z = Lm the 
entrance and exit planes of the amplifying medium, the electric field obeys the 
boundary condition 

F{0,t)^nF{L^,t-At), (6) 
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where TZ is the effective cavity reflectivity ^, At — [L^ — Lin)/wc, and Vc — c/ric 
is the speed of hght within the unloaded part of the cavity, whose refractive 
index is Uc- 

A standard procedure to make the boundary condition isochronous consists 
in defining the new spatial and temporal variables |46| 

Z=^, T = t+^At. (7) 

This transformation amounts to ideally bend the active medium so that its 
entrance and exit coincide, and the delay At accumulated in the trivial prop- 
agation outside the amplifying medium is removed. The boundary conditions 
now read 

F(0,T) =7^F(1,^), (8) 
and the Maxwell-Bloch equations become 

{dT + ^dz)F = KiAP^xF), (9) 

OtP = -f±[FD-{l + id)P] , (10) 
OtD = -f\\[l - D -Re{F*P)] , (11) 

where apsR = 27rc/£c is the cavity free spectral range, 

K=-^(|ln7^2|+a,„L,„), (12) 

is the cavity lincwidth, Cc = n-a\Ln\ + nc{Lc — im) is the optical length of the 
cavity, and 

|ln7^2|+a,,L„, ^ ""^ 

is the adimcnsional pump parameter, while 

< X < 1, is an adimensional loss parameter due to distributed losses inside 
the amplifying medium. In the absence of distributed loss {a^ = 0), x = Oi and 
the pump parameter becomes A = aL^/\ ln7?,'^|, which is its usual definition. 

In the RNGHI for class B lasers it is known that instability occurs when the 
cavity free spectral range is a quantity of order ^/TjjXL (see below). In view 
of this, it is convenient to introduce the new temporal and spatial variables 

r = V^T, C=^Z, (15) 

a 



^TZ accounts for all localized intracavity losses occurring outside the active medium (at 
splices, filters, output couplers, or other components). Namely, li^ = Tif-H2 ■ ■ ■ Ti^, where TZ^ 
is the fraction of power after the j-th lossy element. Distributed losses outside the amplifying 
medium which damp the intensity of the electric field at a rate Oout can also be included in 
the effective reflectivity TZ by multiplying it with Gxp[ — Q^out 
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where 

apsR 27rc , , 

" = , = — : , (16) 

is the scaled free spectral range. If we define also the adimensional decay rates 

, , „ (17) 

the Maxwell-Bloch equations read 

{dr+dc)F = a{AP-xF), (18) 
drP = 7"^ [FD - (1 + iS)P] , (19) 
drD = -f[l - D -Rc{F*P)] , (20) 

with the boundary condition 

FiO,T)^nFiCn.,T), (21) 

where 

Cm = ^ (22) 

a 

is the position of the exit plane of the amplifying medium in the space variable 
C- Note the following useful relation 

2fTCm(l-x) = |ln7^2|. (23) 
2.1.2 Three and four level atoms 

Real lasers are not two-level lasers, even those in which the lasing transition 
can be well described as a two-level system. This is so because pump and relax- 
ation processes connect, incoherently, the two lasing levels to some other levels. 
Although under usual conditions (essentially when the relaxation rates of the 
pumped transition are very large as compared with the rest of relaxation rates) 
this incoherent coupling does not have any influence on the coherent dynamics 
of the system, it turns out to be essential for understanding the meaning of the 
laser parameters, particularly that of the pump parameter A. The connection 
between the two-level theory and the description of three- and four-level lasers 
(which are the usual schemes for describing real lasers) was first treated, to 
the best of our knowledge, by Khanin in his book but he did not analyze 
the consequences of the parameter transformations he derived. Later on, and 
independently of Khanin, these transformations were derived again in |77j ^, 
and their influence on the understanding of laser dynamics was clarified. The 
relation between two-level and three- and four-level lasers is important in our 
context because we will apply our analysis especially to rare-earth doped fibre 

^In lYYI equal relaxation rates for the two lasing levels in four— level lasers were assumed. 
This is not a realistic approximation for NdiYAG lasers. In 1791 . following 1111 . it was shown 
that the same model can be derived even removing that assumption 
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lasers, where the active medium should be properly modelled collec- 
tion of three- or four-level atoms if the dopants are, respectively. Erbium or 
Neodymium atoms ^TE^. 

Eqs. (|18H2U|) . describing two-level lasers, are valid also for three- and four- 
level atoms provided a dependence on the rate of incoherent optical pumping 
Wopt is included in the parameters 7|| and do [771 [7^1 ■ Precisely, we define the 
adimensional pumping rate 

W^^, (24) 

72 

where 72 is the total spontaneous decay rate from the upper level of the lasing 
transition. Then, the parameter transformation necessary for applying Eqs. 
P8H20|I to three- and four-level lasers are 

71! =72(1 + W^), (25) 

and 

W - Sn,3 i^f.. 

where Sn,3 is the Kroenecker S and N is the number of atomic levels. Thus, 
Sn,3 = 1 for three-level atoms, and Si\[,3 = for four-level atoms. 

In fibre lasers the gain parameter a is usually defined in a way different from 
Eq. |(SJ). If (To and Ua are, respectively, the stimulated emission and absorption 
cross sections, one has [7S| 

a = Y [o-o - CTa + (o-o + fTa)do] ■ (27) 

Assuming (Tc = (Ta, the gain coefficient turns out to be proportional to the 
equilibrium total population difference per unit volume JVdo^ as in Eq. (O 

a = (ToA/'do ■ (28) 

Under this approximation, taking into account the definition of the pump pa- 
rameter A, Eq. H13I) . and Eq. (|26|l . the dependence of the pump parameter on 
the pumping rate W can be expressed in the following way 

A = G^-^^\ (29) 

with 

^, Go^JVa^L^. (30) 

Notice that 7|| enters in the normalizations used for writing Eqs. (|18II2U|) . 
and as this parameter is pump dependent in the case of three- and four-level 
lasers, Eq. (|25|l . this must be taken into account when interpreting the results 
derived from Eqs. H18I - I2U|I . The recipe is very simple: (i) replace A by using Eq. 
H29() , and (ii) replace the spatial frequency a by aVl + W in the final expressions 
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(see following sections). This last replacement is due to the normalization of 
the axial coordinate and, consequently, of the spatial frequency, see Eq. lfTB|l . 

The most relevant feature of Eq. H29|) is that the effective two-level pump 
parameter A depends on the actual pump strength in a non-linear fashion 
for both three- and four-level lasers. Although the instability threshold will 
be discussed below, let us comment that for three-level lasers the relationship 
between A and W and the fact that usually G 3> 1 imply that the "factor- 
of-nine"(the ratio of A at the RNGHI instability threshold and at the laser 
threshold) (see Sect. 12.2.21 below) becomes a "factor of (1 + e)"for W, which is 
the actual pump parameter that can be measured in an experiment j77j. This 
is very easy to see: From Eq. (|29|l 

(31) 

and then, the instability to lasing threshold 

W^ins\ _G + 9 G+1 16 

where the limit holds for large G. Contrarily, the case of four-level lasers is very 
similar to that of two-level lasers as in this case 

^0 Al G-9'G-1 G' 

where the limit corresponds to large G again. We see that the instability thresh- 
old for four-level lasers is larger than that of two-level lasers, tending to it for 
very large gain. 

2.1.3 Transverse effects 

So far we have limited our analysis to the plane wave approximation. However, 
transverse effects may play an important role in real experimental situations. In 
fibre lasers, for instance, both the laser and pump field have a transverse spatial 
structure and the doped region has a finite extension [75| . 

To include these elements in the laser equations we first notice that the 
definitions of the dynamical variables F, P and D, Eq. are no longer 

appropriate, because such definitions contain the parameters 7|| and do which 
depend on the pump W. If W is allowed to vary spatially, this introduces an 
undesired spatial dependence in the dynamical variables. The problem can be 
removed by introducing a new definition of F, P and D 

F=-^, P = _2jpi2,/^, D^d. (33) 

The dynamical variables do not depend on W any longer, and they obey the 
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dynamical equations 



{dr + d^)F - a{GP-xF), (34) 
drP = 7"^ [FD - (1 + i5)P] , (35) 
drD = -i[W{l-D)-5N.,3-D-Rc{F*P)] , (36) 

where in a and 7, as well as in the free spectral range a, the decay rate 7|| is 
replaced by 72. 

Now we must include in these equations the dependence of the electric field 
F on the transverse coordinate. We assume that the laser field is a Gaussian 
beam with 1/e radius equal to wq, which is a fairly good approximation j76| . 
This implies that in the Maxwell-Bloch equations all the dynamical variables 
depend also on the radial coordinate r through the scaled coordinate p — r/wg, 
and in the equation for the electric field a term describing diffraction must be 
included 

- V^F +idr+di)F = a (GP- xF) . (37) 

Here V'^ = d^/dp^ + {l/p){d/dp) is the transverse part of the Laplacian and Co 
is the scaled Rayleigh length of the beam i.e. the distance in the propagation 
direction over which the size of the beam varies appreciably. A Gaussian beam 
is described by the function [23 



C) = -TTT exp 

w{0 



+ I — :T-—r I arctan — 



with iiP'iri) = 1 + (C/Co)^5 which is a solution of the empty cavity equation 



(38) 



— VV = acV- (39) 

41,0 

We can write the electric field as F{p, r) — ip{p, C)/(Ci ''') 8.nd, taking into 
account the normalization property of ip 

00 -1 

dppm,p)\'^-, (40) 

valid for any C, we can project the equation for the total field F on the mode 
amplitude / by multiplying each term by ip* and integrating over p from to 
00 

{dr + dc).f = <J (^AG j\pe-p" P-xfY (41) 

Here we have also assumed C *C Co, so that the dependence on C of -i/; can be 
neglected, and ip{p) = exp(— p^). We make the same assumption for the pump 
field, and we denote its maximum amplitude by /S^^^ and its 1/e radius by Wp. 
Introducing the parameter ij = (wo/wp)^, the pump term, which is proportional 
to the intensity of the pump field, can be written as W{p) — P exp{~2rip^). 
Finally, we replace / again with F and introduce the new radial coordinate 
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u — 2(0^. The Maxwell-Bloch equations suitable to describe transverse effects 
read 



{dr + dc_)F = a [g j due-"^^ P - xF] , 
drP = le-^^^FD - {l + iS)P 

drD = ^[(3e-^^{l^D)~5N,3 
with the boundary condition 

f(0,U,T) =7^F(C„„U,T), Cm = 2^/5 



(42) 
(43) 

D-e~''/^Re{F*P)] , (44) 
(45) 



We have fixed the upper hniit in the integral over the radial coordinate equal 
to Um = 2(rni/ifo)^, where is the radius of the amplifying medium. In the 
case of a fibre fm coincides with the radius of the doped region. 

2.1.4 Stationary singlemode solutions 

We now return to the Maxwell-Bloch equations (|18|I - H2U() in the plane wave 
approximation and look for their single frequency (singlemode) solutions, which 
can be found by setting 



F(C,r)=FaC)e-^'^^ P(C,r) = PaC)e- 
The stationary atomic variables are 

1 - iA 



I?(C,t) = A (C) 



1 + A2 + |F, 
1 + A2 



,2 ^ s : 



(46) 
(47) 



Ps = 

where A = 5 — jlu and Fg verifies the differential equation 

dF 
dC 

Using the polar decomposition Fg = Ig^^e''^' and splitting the above equation 
into its real and imaginary parts, one obtains 



di, 

dC 

d0s 
dC 



= 2a 



UJ - 



A 



1 + + F 
aAA 



-X I. 



1 



with the boundary conditions 

/,(0) = 7^2/s(Cm), 

(Cm) - 4>siO) = 2Tm , n integer . 



(48) 
(49) 



(50) 
(51) 
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Integration of Eq. H48|) yields 



[I + In 



hiO A, A-x[l + A2 + J,(C)] 



, . In . 

h (0) X ^ 



= [A-x(l + A2)]2aC. 



(52) 

Particularizing this equation to C = Cm and taking into account Eqs. (|23|l and 
H5U|I . we find that the stationary intensity at the exit of the active medium is 



h (Cm) = ( - 



1 — exp 



1 — TZ'^ exp 



(53) 



It can be easily seen that the lasing threshold, i.e the value of A for which 
Is{C,m) = 0, is 1 + A^. In the hmit of no distributed losses (x ^ 0) Eq. ifS^ 
reduces to 



Is (C„ 



|ln7^2| 



(x = 0) . 



(54) 



The frequency uj of the lasing solution has still to be fixed. To do that we must 
determine how the stationary phase 0s varies along C- If wc insert Eq. H48(l into 
Eq. igni), the latter can be written as 



aAx 



from which it follows that (ps depends on C in a nonlinear way as 



0. (C) - <Ps (0) = (c. - aAx) C - y In III . 



(55) 



(56) 



The total phase shift experienced by the field in a roundtrip from to Cm is 
then 

27r 

'/'s (Cm) ~ <Ps (0) = (t^ - crA) Cm = {uj - (j5 + ajuj) — , (57) 

a 

where we have used Eqs. (|23|l and (|22|) and the fact that A = S — '^uj. Imposing 
the boundary condition (|51|l . wc obtain for the n-th longitudinal mode 



1 



A„ = 



(58) 



7cr 1 + ja 

The singlemode solution associated with the n-th mode has threshold 1 + A^j . 
Recalling that S — (cja — ^^c) /l±, and that jcja — ti-'d < aFSR/2 by definition 
(ttFSR is the cavity free spectral range), we have \S\ < \ict. Hence, the mode 
with the lowest threshold is the one with n = 0, according to our initial choice 
of the reference frequency. For this mode 

Wo — = \ , Ao = :;— = ; . (59) 



1 + 7(7 K + 7_L 1 + 7c + 7_L 

The expression for ljq corresponds to the well-known mode pulling formula. 
From now on A = Ao will be our laser-atoms detuning parameter, and the 
condition of perfectly resonant laser will be equivalent to A = 0. 
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2.1.5 The uniform field limit 

A widely adopted approximation in the study of laser dynamics is the so-called 
Uniform Field Limit (UFL), which is based on the assumption that the steady 
state of the electric field is nearly constant along the propagation direction inside 
the amplifying medium. In other words, in the uniform field limit the electric 
field obeys a periodic boundary condition, as if the cavity mirrors were perfectly 
reflecting. From the analytic point of view this represents a great simplification 
of the problem, because it allows for an expansion of the electric field in terms 
of the Fourier modes exp(zQ!„C) of the empty cavity, where an = na is the 
frequency of the n-th mode and the singlemode stationary solution associated 
with mode tt, = is spatially homogeneous. From the numerical point of view, 
the benefits of the UFL are even more relevant, because the expansion of electric 
field in Fourier modes allows to transform the original set of Maxwell-Bloch 
equations from partial differential equations to ordinary differential equations, 
which can be numerically integrated much easier and faster. Although our work 
on the RNGHI shows that the above mentioned features of the UFL persist 
even outside that limit for class B lasers (see Sect. I2.3|) . in the most complex 
situations - for instance when transverse effects or inhomogeneous broadening 
are included in the model - we still limit our analysis to the UFL for the sake 
of algebraic simplicity. 

In |8(Jj we introduced a new technique for dealing with the UFL. In compar- 
ison to the standard technique [SSI ' technique is advantageous because 
it allows to clearly demonstrate that the only requirement for applying the UFL 
is that the resonator reflectivity must be close enough to unity (the single-pass 
gain does not have to verify any constraint). In (SOI the simplest case of exact 
resonance and absence of distributed loss was considered. In this section we 
show, following [SOI, how the Maxwell-Bloch equations in the UFL can be de- 
rived for homogeneously broadened lasers in the plane wave approximation in 
the presence of both detuning and distributed loss. The generalization to other 
situations (including, e.g., transverse effects or inhomogeneous broadening) is 
straightforward . 

To derive the dynamical equations in the UFL we first observe that in the 
limit TZ^ 1 the solution of Eq. H52|l for the n-th mode varies slowly along ( 
according to the equation 



where T = 1 — TZ^ is a small parameter. We thus see that, independently of the 
value of all other laser parameters, the laser intensity is almost uniform along 
the active medium in the limit TZ^ —> 1, hence the name UFL. If we insert the 
above expression for Is{C) in Eq- fSH|l . and take into account Eqs. and 
H58|l . we find that the stationary phase also varies linearly along ( according to 
the simple expression 




(60) 




(61) 
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which shows that the space frequency associated with the n-th solution is na. 
We focus on the n — solution for which the stationary phase is constant along 
C,. Thus, in the UFL, this stationary solution can be written as 

pUFL _ 1-iA pUFL nUFL _ l+A'' l"^/ 

— A ; — A ■ 

Now we introduce new dynamical variables obtained from the old ones dividing 
them by the stationary solution n = outside the UFL and multiplying by the 
same solution in the UFL 

F'iCr) = ^^§^F^-\ (63) 

^^^^ = TTa^TXIo = i + A^ + /s(C) • ^''^ 

Since i^s(O) = TZFs{Cm), the new electric field F'{(, t) obeys a periodic boundary 
condition 

F'{0,T)^F'iCm,r). (67) 
The Maxwell-Bloch equations 118I - I20|I for the primed variables take the form 
{dr+dc)F' = aT]iC)[AP' -{l-iA)F'] , (68) 
drP' = 7"^ [F'D' -{1 + iA) P'] , (69) 



with 



drD' = 7 



(70) 



Because of the periodicity condition H67|) these equations are particularly suit- 
able for numerical integration, as we shall discuss in Sect. 12.3.51 
Obviously, Eqs. H68II70|I admit as a stationary solution the UFL stationary solu- 
tion H62|l . Yet, they are still equivalent to the original Maxwell-Bloch equations, 
and this fact is reflected in the dependence on C of and I^. But from Eqs. 
P7|) . and ^ we see that in the UFL 

^^(0 - ^ ' t(^-I-) +0{T') . (71) 



4"^^ VCm 2 

viC) = i + ^''/^V (i--0+0(r^), (72) 

Thus, in the very limit we can set both terms equal to unity, and, dropping the 
primes, the laser equations in the UFL take the well-known form 

{dr + dc)F = a[AP-{l-iA)F] , (73) 

drP = j-'^[FD-{l + iA)P] , (74) 

drD = "/[I- D-Re{F*P)] , (75) 
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with the boundary condition 



F(0,t) =F(C„,t) , 



.m — 



27r/a. 



(76) 



These equations differ from the general ones (Eqs. (|18|I - H21() ') only in one re- 
spect: The (localized) cavity losses, represented by the mirror reflectivity 7?,, 
have disappeared from the boundary condition (which now is periodic) and 
have appeared in the dynamical equation for the electric field through the term 
—aF, which contains both localized and distributed losses. The appearance of a 
detuning in the equation for F as well as the replacement of 6 with A in the 
equation for P are only apparent changes, due to the fact that in Eqs. H73l - I75|l 
the reference frequency is the lasing frequency luo, while in Eqs. H18|) - (|20|l the 
reference frequency is the empty cavity frequency tOc- 

Let us stress that to derive the Maxwell-Bloch equations in the UFL we 
had only to assume that TZ^ 1. Eqs. (|71H72|I show that this assumption 
suffices to approximate Is{C) / 1^^^ and 7]{C) with unity, for any value of the 
pump parameter A greater than one, no matter how large it is. Therefore, 
the condition of small single pass gain aL^ — > 0, which invariably accompanies 
the condition TZ^ — > 1 in the literature about the UFL, is actually completely 
superfluous. Notice also that the amount of distributed loss does not influence 
the vahdity of the UFL. 



2.2 The RNGH instability in the uniform field limit 



The uniform field limit (UFL), rigorously introduced in Sect. 12.1.51 is the case 
in which the RNGHI can be analyzed in the simplest way, as well as the limit 
where the effects of additional features not considered in the original analyzes by 
Risken and Nummedal ^22^ and by Graham and Haken can be studied most 
easily. Still another reason for the interest in the UFL is that it still captures 
the very basic signatures of the multimode instability of ring lasers found in 
more complex models. 

In this section, we shall first derive again the classical results of RNGH (Sect. 
I2.2.1|l . Then we shall consider the relevant case when the population inversion 
is an extremely slow variable (class B laser). This applies, in particular, to fiber 
lasers (Sees. 12.2.21 and 12.2.31) . Next we will show that detuning plays no role 
in class B lasers (Sect. 12.2.4(1 . Finally, we shall analyze the spatial (transverse) 
effects due to the modal structure of fiber lasers (Sect. I2.2.5|l . 

2.2.1 The RNGH instability 

The original studies by RNGH considered a perfectly resonant two-level 

ring laser in the UFL. The Maxwell-Bloch equations that describe such laser 
are given by Eqs. introduced in Sect. 12. 1.51 with A = 



{dr + dc)F 

drP 
drD 



g{AP - F) , 
^-\FD~P), 
^[l-D-Rc{F*P)] , 



(77) 
(78) 
(79) 
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supplemented by the periodic boundary condition 

F(0,t) = F(C„„t) , Cm = 2^5. (80) 
The steady, spatially uniform, lasing solution is given by Eq. I(j2|l with A = 



Ps - 



El 

A 



A' 



(81) 



where 4> is an arbitrary phase, and it exists for ^ > 1 (A = 1 corresponds to 
the first laser threshold). This solution is the basic, singlemode lasing solution. 
Other singlemode solutions exist (in fact an infinite countable set) that can 
be obtained from Eqs. (f77jl - l|5n|l by setting F {(, t) = Fg.n exp [i {na(^ — a;„T)], 
P (C, t) = Ps,n exp [i {na( — w„r)], D {(, r) = Ds^n. These solutions were com- 
puted in Sect. 12.1.51 outside the UFL and correspond to different values of the 
longitudinal index n. As we showed in that section, all these additional single- 
mode solutions have a larger oscillation threshold and hence solution Eq. H81|l 
is the first lasing mode. 

Stability of this singlemode solution is analyzed by perturbing the state of 
the system as 



F(C,r) = F + e'-^SFiCr)., 
P(C,r) = P, + e^*<5P(C,T)^ 
DiC,T) = D, + SD{C,t), 



(82) 
(83) 
(84) 



and linearizing the dynamical equations for the perturbations. The study is 
facilitated by expressing any perturbation SX in terms of its real (SX^) and 
imaginary (dXi) parts — note that 5D is real by definition. The system of 
equations can be expressed in vector form as 



drSX ^L-SX, 
where 6X = col (JFr, (5Pr, 5D, 5Fi, 5Pi) and 



L 



(85) 



— <J — 9(; 


aA 























-l\Ps\ 


-l\Fs\ 


-7 

















—<7 — 


aA 
















(86) 



with Fs, Ps, and Ds given by Eq. l|HT|l . 

Because of the linear nature of the system (|85|) . any solution 5X can be 
calculated as 5X {(, t) = J2\ SX\ {() c^'^ . On the other hand, given the depen- 
dence of L on the spatial coordinate C only through the gradient 5^ , any SX\ (C) 
can be written as SX\ = J^a ^^>;a exp {ia(). (Outside the UFL, this ansatz 
is not valid; see below.) Hence any solution 6X can be finally calculated as 



SX{C,t) 



SX 



A,ae 



17 



Thus the dynamical problem (|85|l is transformed into the following eigenvalue 
problem 

where La is given by matrix L with substituted by ia. Given the block- 
diagonal form of La the characteristic polynomial V (A; a) that determines the 
Lyapunov exponents A is factorized as 

P (A; a) = Vi (A; a) Vr (A; a) , (87) 

Pi (A; a) = + (j^^ + a + ia) X + i-f-'^a , (88) 
VR{X;a) ^ X^ + {j-^+j + a + ia)X^ + 

[A + -fa + ia (7"^ + 7)] A + 2a {A - 1) + iAa . (89) 

These polynomials, equated to zero, allow to determine the dependence A (a) of 
the eigenvalues on the spatial frequency a of the sidemode perturbation. The 
boundaries of the unstable domain for the multimode instability correspond 
to cases where ReA (a ^ 0) = (for a — one is considering the well known 
singlemode, or Lorenz-Haken, laser instability, which is not treated here). The 
singlemode solution will become unstable against sidemodes of spatial frequency 
a if ReA (a) > 0. 

It is easy to show that Vi has no roots with ReA (a) > for any a. Only 
marginally A (a = 0) =0 is always a root. This solution does not entail an 
instability as it never gets positive real part; in fact that root is merely reflecting 
the phase arbitrariness of the lasing solution. As Vi is associated with the 
subspace formed by col {SFi,S Pi), see Eqs. (|S5|) and (|SS1 ). and that subspace 
controls the possible growth of perturbations in phase-quadrature with the lasing 
mode, one concludes that the RNGHI is not a phase instability but an amplitude 
instability l^nHSOl- This result remains valid outside the uniform field limit, as 
well as when other factors (inhomogeneous broadening, etc) are included, as far 
as the resonance condition is maintained. As we discuss below (Sect. 2.2.4), 
detuning does not have any influence in the stability of class-B lasers, which 
are our main interest. Nevertheless, in class-A and class-C lasers, detuning has 
a large influence on the stability properties |41l 1431 ICT] . 

Regarding the polynomial Vr — which, mutatis mutandi, governs amplitude 
instabilities — the boundaries of the multimode instability are found by setting 
A = —iu). Upon splitting the thus obtained polynomial (now in w) into its real 
and imaginary parts and by equating them to zero one can solve for a and lo as 

a± = u;±(l + ^^^], (90) 



ljI = ^ 3{A-1)--/^±VR , (91) 
R = {A-l)iA-9)-6j^iA-l)+-f\ (92) 

The homogeneous solution turns out to be unstable for values of a verifying 
\ct-\ l£ \(^\ 5: This unstable domain has the shape of a tongue in the plane 
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Figure 1: RNGHI threshold in the (r, a) plane for the three values of 7 marked 
in the figure. 

< A,a >, and the two branches a± merge at the critical point {Ac, etc) defined 
by the condition R = 0. The critical pump represents the minimum value of A 
for which the instability can exist, and is given by [^1^ 



The value of ac is obtained from Eqs. H9U|I and H91|) by setting R = and 
A = Ac- The results of this linear stability analysis are summarized in Fig. ^ 

2.2.2 The RNGH instability in class B lasers 

Class B lasers are defined by the inequalities 7|| <C k <C 7_l concerning the decay 
rates of the population inversion, intracavity field amplitude, and medium polar- 
ization, respectively. With the used normalizations those inequalities transform 
into 7 <C cr ^ 7~^. This limit is especially interesting as fiber lasers — and 
erbium-doped fiber lasers in particular — belong to this class. On the other 
hand that limit is interesting also from the theoretical viewpoint as it allows 
to obtain a wealth of analytical information. The latter does not show up in 
the "classical" version of the RNGHI we have just considered, as the linear sta- 
bility analysis of the singlemode solution against multimode perturbations is 
completely analytical. In this section we apply the class B limit to the general 
expressions already obtained above. In the following section we show how the 
same results can be derived using asymptotic techniques that employ 7 as an 
expansion parameter. 




(93) 



19 



The expressions for the instabihty boundaries in the class B hmit are re- 
trieved from Eqs. ipn|) ~ (|^ by setting 7 = 



al = \\i(A-l)± v/(A - 1) (A - 9)J , (94) 
L0± = a±. (95) 

The critical pump is obtained from Eq. H93|) by making 7 = and reads 

A, = 9, (96) 

which is in fact the minimum value A^ can attain. As the first (lasing) threshold 
occurs at a pump A= 1, we conclude that in order to have a RNGHI the pump 
parameter A should be at least a factor of nine greater than its value at lasing 
threshold. This is the famous "factor-of-nine" . 

As for the sidemode critical frequency ac, we obtain 

ac = %/T2. (97) 

So far we have assumed implicitly a continuum of longitudinal modes la- 
belled by their wavenumber offset a. But it must be remembered that the 
periodic boundary condition l(5n|l imposes, in particular, that the perturbation 
wavenumber a must be an integer multiple of the scaled cavity free spectral 
range a, Eq. l(T?)|l 

2-KC . 

a = Un = na = n , n mteger. (98) 

Thus, in order to have an instability at the lowest possible pump A — A^ the 
critical frequency ac must verify H98|l : alternatively, the cavity (optical) length 
Cc must verify 

TTC 

Cc = n — , n integer. (99) 

\/37||7-L 

Hence the shortest cavity length that allows the RNGHI a,t A = Ac is 

Cc,min = ^ , (100) 

which is the "critical cavity length"which we referred to in the introduction. 

Next we analyze the infiuence of the actual value of the cavity length (or, 
alternatively, the consequences on the instability of discrete nature of the cavity 
modes) on the instability. The analysis will clarify the meaning of £c,min ) which 
we anticipate is not a true minimum (critical) value for the cavity length, but 
only the minimum value of Cc for the instability to occur at the lowest pump. 
Let us consider the RNGHI boundary, Eq. (|94|1 . shown in Fig. |21 Besides 
the critical point {Ac,ac) another relevant point is the minimum of a_ vs. A, 
which is easily determined from Eq. 194() and evaluates to Omin = 2-1- \/2- This 
quantity is relevant as it settles the minimum frequency of a sideband that can 
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Figure 2: Enlargement of Fig. 1 for 7 = 0. 

be linearly amplified after the RNGHI. If the cavity is short in the sense that 
its free spectral range a is greater than amin (i.e. if Cc < Ci, = — =£=) 
then the sideband with lowest threshold (which we denote by Athi) will be the 
first one (the one with a = a). This threshold value is obtained from Eq. H94|l 
and reads 

Athr = ^ (2 + _ ^4 _ i2a2 + a*) , (101) 

where a must be substituted by a and the expression is only valid for a > 
amin = 2 + ^. As by decreasing Cc, a increases, the RNGHI is produced at 
increasing pump values according to Eq. (|101|l so that Athr — * 00 for Cc — > 0. 
This reasoning evidences that £c,min is not a true minimum value for the cavity 
length in order to observe the RNGHI, contrarily to what is commonly believed. 
Clearly, by lowering Cc Athr can become so huge that, in practice, the RNGHI 
can be ruled out. But we stress that this is a practical, not a fundamental, 
limitation. On the other hand when Cc > Ci, then a < amin and the instability 
is produced at a higher order sideband. The identification of the sideband with 
lowest threshold is now a more involved task (see ;77^ for a discussion). In 
general, for each value of Cc, Athr must be computed from Eq. (|101|l by taking 
a = na with n = 1,2,3,... and the lowest value of A thr corresponds to the 
actual instability threshold value. Fig. 13 displays the actual RNGHI threshold 
as a function of the cavity length. We note that the graph has periodically 
recurring minima. This is due to the fact that whenever Eq. (|99|l is verified, 
Athr = Ac = 9. All the above analysis can be applied to three- and four- 

^We note that this equation corrects some typos appearing in Eq. (45) of I77| . 
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Figure 3: RNGHI threshold as a function of the cavity length for 7 = 0. 



level lasers making use of the transformations defined in Sect. 12.1.21 (i) replace 
A by using Eq. (|29|l . and (ii) replace the spatial frequency a by a\/l + W 
in the final expressions (remember that W is the actual pumping strength). 
As already advanced in Sect. 12.1.21 for four-level lasers the situation is quite 
similar to that of two-level lasers, but for three-level lasers the scenario changes 
dramatically as we discuss next. On the one hand the critical value for two- 
level lasers etc = \/T2, Eq. l|!?7l) . must be replaced by y^l2 (1 + Wc) where Wc 
is the minimum (critical) instability threshold referred to the actual optical 
pumping W, which is given by 

G + 9 

Here G is the gain coefficient for three-level lasers defined in Eq. H30II . For 
G 1, as is typical in erbium-doped fiber lasers fTH, Wc 1 and then ctc^'^'' = 
1/12 (1 + Wc) ^24, which is a factor larger than for two-level lasers. 
With respect to £c,min 



risL) TTc G-9 nc 

^6711 7_L V G V^7||7-i 

This value of £c,min is times smaller than the corresponding two-level value 
H1U(J|) . However the most fascinating difference with two-level lasers is the ex- 
tremely low instability threshold obtained in three-level lasers. Denoting by Won 
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the lasing threshold referred to W, one has Wc/Won — 1 + 16/G when G ^ 1, 
as already discussed in Sect. 12.1.21 fsee Eq. Thus in erbium-doped fiber 

lasers the RNGHI is predicted to occur just above threshold |75[ I77| . Finally, 
the variation of the instability threshold with the cavity length in three-level 
lasers is much less pronounced than in two-level lasers and, in fact, the cav- 
ity can become substantially shorter than £c,min while the instability threshold 
keeps moderate values; see Fig. 6 in |77] . 

2.2.3 Asymptotic expansions for class B lasers 

The results discussed in the previous section have been obtained by applying the 
class B limit 7 to the general expressions obtained in Sect. 12.2.11 which are 
valid for arbitrary values of the laser parameters. In other instances (e.g., outside 
the uniform field limit) this strategy is not possible as general expressions are not 
available. In such cases one can still obtain a lot of analytical information if the 
study is done, ab initio, for class B lasers. Those treatments rely on asymptotic 
expansions of the problem that use 7 as a smallness parameter. Although in the 
UFL this analysis is not necessary, we prefer to introduce the technique at this 
point as here the explanations are more transparent and straightforward than 
in more complex cases that will be considered below. 

The starting point of the analysis is the characteristic polynomial governing 
the RNGHI, which must be equated to zero in order to obtain the eigenvalues 
A that govern the stability of the singlemode lasing solution. In the UFL it is 
given by Eq. H89() . which we recall for convenience 

VR{X;a) = X^ + {j-^+j + (T + ia)X^ + 

[A + 7(7 + ia (7"^ + 7)] A + 2(7 (A - 1) + iAa . 

(We ignore the polynomial Vi associated with the phase as we showed it does 
not contain any instability.) Now, as < 7 ^ 1 defines class B lasers (together 
with (7 ^ 7"), we assume the following ansatz 

A ^ Ao +7A1 +7^-^2 + ■ • ■ , (104) 

for the eigenvalues. The ansatz is substituted into Vb. (A; a) and the resulting 
polynomial is expanded in series of 7 as 

00 

n=N 

and, in our case, the expansion ljlQ5(l starts at N = —1. 

Now Vn must be equated to zero. As the expansion is assumed to be uni- 
formly valid for any 7, we impose that each of the p„ be null. Starting at the 
leading order n = iV (= — 1), we have 

p-i = Ao (Ao -I- ia) , 
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which has the roots 



and A| 



'0 



(106) 



The next order [n 



0) reads 



Po = AoAi + (Ao + ia) {Xi+A + A^) + 2a {A ~ 1) + crAj]. 



(107) 



Making use of the first root Aq = and equating po to zero we obtain Ai = 
—A + 2ia {A — 1) /a. As ReAi < the roots corresponding to Aq = do not 
entail an instabihty and can be discarded. Then we must seek instabihties 
associated with Aq = —ia, Eq. H1U6|) . Substituting this root into Eq. (|107f) and 
making pq = 0, we obtain 



As ReAi = we must continue the analysis. At the next order (n = 1), once 
Eqs. (|106|1 and (|108|l have been used, and after making pi = we obtain 



The analysis can stop here as ReA2 is not identically zero. The instability 
boundary ReA = reads in this case ReA2 = 0, which exactly reduces to the 
instability boundary for class B lasers Eq. H94I) . or Eq. (|101|l . as can be checked 
easily. Regarding the oscillation frequency at the instability threshold, which 
we denoted by uj in the previous sections, it corresponds to ImA. To the leading 
order ImA = ImAo = —a, Eq. H106() . Thus it is predicted that lu = a, in. 
agreement with Eq. H95|l . 

This kind of analysis will be used later in order to treat more involved 
problems which, usually, do not admit general analytical expressions to which 
the class B limit 7^0 can be applied. 

To conclude, we note that the above expansion assumes implicitly that all 
parameters are of order 7'^. Obviously, additional scalings can be incorporated 
if needed, as it occurs in the next section. For instance, in order to fully under- 
stand the bifurcation in the presence of inhomogeneous broadening, the scaling 
a = 7^^a_i must be studied separately jSlj. 

2.2.4 Role of detuning in class B lasers 

Throughout this paper we shall be dealing with resonant models, as we have 
just done in the previous sections. Here we prove that detuning between the 
cavity and the gain line, which is almost unavoidable in real experiments, can 
be ignored in class B lasers. This is important as it demonstrates that the 
predictions of resonant models remain valid, to the leading order, even when 
cavity detuning is present. The analysis will be done in the UFL for the sake of 



Ai = iaa 1 ^ {A - 1) 



(108) 



a- 




(109) 
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simplicity. Thus the starting point of the analysis are Eqs. (|73|I -H75 (I 



{dr + d<:)F = a[AP-{l-iA)F] , 

drP = -f-'^[FD^{l + iA)P] , 
drD = 7 [1 - L> - Rc(F*P)] , 

with the boundary condition H7t)|) 

FiO,T)=F{U,T) , U^2tt/&. 

where a is the scaled cavity free spectral range, Eq. l(TB|) . and 

A = -^. (110) 
1 + 7Cr 

As pointed out in the derivation of the equations in the UFL, by definition 
\S\ < ja/2. In the class B limit 7 < 1 and a ^ O (7''), hence |A| < ^ja to the 
leading order. Next we determine the maximum admissible order of magnitude 
of A. For that we note, based on the results of the resonant model discussed 
in the previous sections, that a must be at most a quantity of order 7" for the 
instability to occur at physically realizable pumping levels. Hence A is, at most, 
of order 7. This is shown by the following argument: As the case of small a 
(say a ^ -f^ with k — 1,2,...) poses no problem (in this case A would be still 
smaller) we concentrate on the large free spectral range limit a 3> 1. In this 
case, as discussed in Sect. 12.2.21 the first sideband (a — a) is the one with 
lowest instability threshold, whose value is given by Eq. H101|) . The asymptotic 
form of such equation in the limit a 3> 1 reads 

Athi- > 2^ ■ 

Then, if a ~ 7^^, Athr ^ 7^^ ^ 10^° (for typical fiber lasers 7 ^ lO^'^), and 
this is nonsense^. The conclusion is that, unless completely unrealistic pumping 
values are considered, a must be a quantity of order 7° (at most) and then an 
appropriate scaling for class B lasers is A = 7A1, with Ai of order 7°. 

The singlemode solution with lowest threshold is the stationary solution 
given by Eq. Ht)2|) . that we recall here for convenience 

F, ,^ . . „ 1 + A2 



F, = ^A-l-A^e'^ , Ps = (1 - »A) , D 



A ' ' ' A 

For A = this solution reduces to the resonant singlemode solution (|81|l ana- 
lyzed in the previous sections. The lasing threshold 1 + A^ differ from that of 
a perfectly resonant laser only by a term of order 7^ as A is of order 7. 

The stability of the stationary solution can be studied as in Sect. 12.2.11 
In this case however, the presence of the atomic detuning A introduces an 

^Even for a moderate value as a = 20 (which is still of order 7*^), the actual instability 
threshold is Ajij^ = 202.01, which is clearly a huge pumping ratio. 
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entanglement between the real and imaginary parts of the fluctuations, and 
the fifth order characteristic polynomial cannot be any longer factorized as in 
the resonant case, Eq. H87|) . The analysis of the eigenvalues is much more 
complicated and it is convenient to apply the asymptotic expansion presented in 
Sect. 12.2.31 We do not give the details of the derivation, which is straightforward 
but lengthy. We just want to emphasize that the relevant eigenvalue obtained 
in this way reads A = —ia + 7A1 + 7^A2 + O (7'^) , where Ai and A2 are the very 
same we obtained in the resonant case analyzed in Sect. 12.2.31 In particular this 
means that the detuning A = 7A1 does not have any influence to the leading 
order (probably it appears at corrections of order 7^ or smaller, but these do 
not control the instability). The conclusion is then clear: In class B lasers the 
influence of the cavity detuning does not show up to the leading order. Then 
one can safely ignore that detuning and the results obtained in resonant models 
can be considered as an excellent approximation to detuned models. 

2.2.5 Spatial efTects 

When the plane-wave approximation is abandoned and one assumes, more re- 
alistically, that the laser operates on the fundamental Gaussian mode (TEMoo), 
the RNGHI (as well as the Lorenz-Hakcn instability) turns out to be strongly 
influenced by the relative size of the beam and the amplifying medium. If rm 
is the radius of the amplifying medium and wq the beam waist, the relevant 
parameter is Um = 2(rni/wo)'^. In two-level lasers it was demonstrated that the 
instability disappears in the limit Um ^ 1 [12 EH |^ , but it is recovered as 
soon as Wm becomes of order unity [SSJ . This result can be easily under- 
stood considering that in the limit Um <C 1 of very narrow amplifying medium 
the model equations for a Gaussian beam reduce to those of the plane-wave 
approximation, because the wavefront can be assumed to be plane inside the 
amplifying medium. The results of |86| apply well to fibre lasers, where the 
dopant is usually confined inside the fibre core, whose radius rcoro is typically 
smaller than the waist of the laser beam. Yet, we know that when considering 
fibre lasers one must distinguish between three- and four-level atoms, because 
of the different dependence on the pump of the equilibrium population difference 
do- 

The Maxwell-Bloch equations for a Gaussian beam in three- and four-level 
lasers are Eqs. (|42H44|I . derived in Sect. 12.1.31 Here we limit our analysis to the 
UFL. Following the same steps as in Sect. 12.1.51 it can be easily shown that the 
Maxwell-Bloch equations in the resonant case in the UFL are |79l W!\ 




(111) 



(112) 



7 /3e"''"(l-i?)-5^,3-^-e 



(113) 



26 



8 



6 - 



4 - 



2 - 



0.0 




0.2 



0.4 



0.6 



0.8 



1.0 



X 



dop 



Figure 4: Lasing threshold j3i and minimum instability threshold P2,min as 
functions of the scaled radius of the doped region Xdop for three values of the 
gain parameter G: G = 600 (solid lines), G = 200 (dashed-dotted lines) and 
G = 60 (dotted lines). For each value of G the singlemode solution exists and 
is stable in the region between the two lines. 



with the periodic boundary condition 



F{Q,u,t)^F{C„„u,t), Cm = 2tt/& . 



(114) 



Here we assumed without loss of generality that both F and P are real. We 
recall that /3 is the intensity of the pump field, and the parameter rj is defined as 
rj = {wq/w-pY, where w-p is the waist of the pump beam. The parameter 5n,3, 
which is equal to 1 for three-level lasers and to for four-level lasers, makes 
the difference between the two situations. 

Transverse effects in Erbium-doped (three-level) fibre lasers were investi- 
gated in 'K7'. It was assumed that the dopant is confined inside the core of the 
fibre, i.e. < rcorc- The ratio 



•^dop 



(115) 



was used as a free parameter, ranging to 1. For fibre lasers of the type used 
in |72| adequate parameters are Wp = 2.086 ^m and wq — 2.828 /im, which 
impHes rj = 1.839, and rcorc — 1.875 /xm [721121]. Since u^i can be written 
as Um ~ 2(rcorc/wo)^a;^Qp, with our parameters we have Um — 0.879 aj^^p. 
The proper pump parameter now is the pump intensity /3. We call Pi the 
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Figure 5: Solid lines: output intensity as a function of the scaled radius of the 
doped region Xdop for G = 200 and for different values of the pump intensity, 
(3 = 1.25, 1.5, 1.75, 2.5 from bottom to top. Dashed line: output intensity in 
correspondence of the minimum instability threshold /32,min- 



lasing threshold and (32 the threshold for the RNGHI. The minimum instability 
threshold, which is the one corresponding to the critical frequency space Uc 
where the upper and lower instability boundaries merge, is denoted by /32.min- 

A first important result is shown in Fig. ^ where the dependence of Pi and 
/32.min On the normalized doped radius x^op is shown for three experimentally 
accessible values of the gain parameter G . All curves display a vertical asymp- 
tote for small doped regions, and it may be noted that, for each value of G, the 
asymptote for /32,min corresponds to a value of x^op which is about three times 
larger than the one associated with the asymptote for (3i . The existence of these 
asymptotes can be explained making reference to the plane-wave limit, to which 
the model reduces in the limit Xdop — *■ (which implies Um —> 0) |87| . It can be 
shown that the two asymptotes are well approximated by a;dop = l/vCK879G 
and Xdop — 3/\/0.879G, where the factor 0.879 comes from Um = 0.879 a;jQp. 

Thus, there exists a band of values of Xdop between the two asymptotes, 
where singlemode lasing is stable for all pump strengths. One is tempted to 
conclude that to stabilize the laser it suffices to dope only a small fraction of 
the fibre core. But Fig. ^ also tells us that the lasing threshold fix rapidly 
increases for small Xdop, which means that, for a given pump the output 
intensity decreases. A compromise between stability and not too low output 
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intensity must be found. 

Fig. [S] shows how the output intensity / depends on a;dop for four increas- 
ing values of the pump intensity (3, and G — 200. All curves are similar and 
display a maximum that shifts to larger doped area radii as (3 increases. Also 
displayed, as a dashed line, is the critical intensity /2,min associated with the 
minimum instability threshold /32,min- All stationary intensities in Fig. [Slwhich 
are above /2,min are unstable (multimode emission takes place). We show only 
the case G — 200 for the sake of clarity. The vertical asymptote of /2,min cor- 
responds to that of /32,min in Fig. ^for G — 200. Interestingly, /2,min shows a 
relatively large plateau for 0.4 < Xdop ^ 0.7. In this region the maximum stable 
singlemode intensity available is almost insensitive to Xdop, what constitutes a 
genuine transverse effect, not explainable in the plane- wave limit. 

The best operating conditions are given by an interplay between plane- 
wave limit and transverse effects. In the particular case considered (G = 200), 
the plane-wave limit indicates that stable singlemode emission is guaranteed if 
0.1 < 2^dop ^ 0.25, but in this region very large values of the pump intensity 
are required in order to obtain a large output intensity. However, the particu- 
lar dependence of the output intensity on a;dop, which is a consequence of the 
transverse profile of the laser and pump beams, shows that stable emission with 
sufficiently high intensity can be achieved for smaller pump intensities when the 
radius of the doped area is larger (xdop — 0.6). 

A similar analysis was performed in Ref. |^ for a ring Nd:YAG (four-level) 
laser using Eqs. (|llimi3|l with 5n,z = 0. At difference with the erbium-doped 
fiber laser, in the case of a bulk Nd:YAG laser the dopant distribution can be 
assumed to be homogeneous, and to some approximation the same holds for 
the pump (at least in side-pumped designs). Then, the only transverse effects 
that have to be taken into account are the Gaussian profile of the laser field 
and the finite dimension of the active medium (the laser rod). In |7H1 it was 
demonstrated that the instability occurs whenever the beam waist wq is slightly 
larger, or at least equal to r,„, the radius of the active medium's rod. To 
accomplish this is not a trivial task, as it is discussed in detail in [7^1 , but it 
seems feasible with present day technology. As Nd is a four-level medium, in 
this case the instability to lasing threshold ratio is slightly larger than 9 in the 
optimal conditions (see Sect. 12.1.21) . i.e., in this case the instability threshold is 
well separated from the lasing one. We refer the reader to '79' for full details. 

2.3 Outside the uniform field limit 

The study of lasers instabilities outside the uniform field limit is in general a 
difficult task, because even the stationary state is spatially dependent jSHl- But 
the best candidates for the observation of the RNGHI are fibre lasers, for which 
the UFL usually does not apply, because the effective refiectivity parameter TZ 
is far from unity. On the other hand, these lasers satisfy very well the class-B 
limit, because the parameter 7 is as small as 10~^ .22, .22, .89,- This is a 
fortunate circumstance because, as we shall show in the present section, at least 
in the study of the RNGHI, the class-B limit is as powerful as the UFL, in 
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the sense that the analytical results obtained in the UFL can be generalized to 
arbitrary TZ in class-B lasers. 

Specifically, in this section we will show that for class-B lasers outside the 
UFL the following results hold: (i) the boundaries of the stability domain of 
the homogeneous solution can be determined analytically, (ii) the simple self- 
pulsing solutions can be studied in a completely analytic way, (iii) the multimode 
dynamics arising from the instability can be described in terms of a limited 
number of Fourier modes. In all cases cavity detuning will be ignored after the 
discussion in Sect. 12.2.41 

This section is organized as follows: In Sect. I2.3.1l we derive the RNGHI out- 
side the UFL. Then we generalize the result by taking into account distributed 
losses (Sect. KT^ . In Sect. we derive generalized rate equations for 

class-B lasers, and in Sect. 12.3. 4l we use those equations to analyze self-pulsing. 
In Sect. 12.3.51 we discuss the issue of the numerical integration of the model 
and present some numerical results in order to compare the different integra- 
tion methods. Finally, in Sect. 12.3.61 we address the issue of the super- or 
sub-criticality of the bifurcation. 



2.3.1 RNGH outside the UFL 

First we consider the case when no distributed losses exist in the laser cavity 
(x = 0) 8.nd recall the Maxwell-Bloch equations H18I) ~ (|2U|I valid in the resonant 
case 

{dr + dc)F = (jAP, (116) 
drP = 7"^ {FD ~ P) , (117) 
drD = -f{l-D-FP), (118) 

together with the boundary condition 

F(0,r) =7^F(C„„r), Cm = 27r/<5 , (119) 

where C = and ( — Cm denote the location of the entry and exit faces of the 
active medium, and a is the scaled cavity free spectral range, Eq. H16|l . 

In these equations we have assumed that the fields F and P are real because 
we showed in Sect. I2.2.1l that in the resonant case and in the UFL, the RNGHI 
is an amplitude instability. Lifting the UFL does not change the nature of the 
bifurcation, as can be checked following the lines of our study in Sect. 12.2.11 

Unlike in the UFL, now all variables at steady state are in general functions 
of the coordinate C, as we showed in Section r2.1.4l . According to Eq. H54|l the 
steady value of the field intensity /g = | | at the medium exit face is given by 

/.(Cm)-^^(A-l), (120) 

and A — 1 corresponds to the first laser threshold. 

As in previous sections, the linear stability analysis is performed by writing 
F (C, t) ~ Fs (C) + SF (C) exp (Ar), and similar expressions for P and D, and 
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linearizing the resulting dynamical equations. As now the problem depends 
explicitly on the longitudinal coordinate ^ and the boundary condition is not 
periodic, one cannot make the ansatz 6F ^ exp (laC) ^.s in the UFL. In fact, now 
the spatial dependence of 6F is unknown a priori and an equation for it must be 
sought. After solving for the material perturbations we set SF{() = Fs{()Sf{() 
and obtain the following equation for Sf 



dC 

H{Is) 

whose formal integration yields 



-\Sf-H{QSf, 

a A 7A(A + 7) + 
1 + /s (7 + A) (1 + 7A) + jls ' 



SfiO = 5/(0)exphAC-V^(C)] , 



dCH {Is iO) 



The evaluation of the integral tp {() requires a change of variable from C to Is 
through the steady state equation dIs/dC, = 2aAIs/ {\ + Is), after which 
the following result is obtained [TO] 



7A 



■In 



^s(C) , 2 + 7A ^^ J1 + 7A) (A + 7)+7/s(C) ^ (^21) 



2(1 + 7A) 7,(0) 2(1+7A) (1 + 7A) (A + 7) + 7/, (0) 
Finally we return to the initial perturbation 5F{C^) and write 

F (C) 

6F (C) = 5F (0) exp [-AC - (C)] • 



(122) 



The boundary condition H119|) applied to the perturbation SF (C) gives the fol- 
lowing characteristic equation 



A 



-ma - 



V' (Cm) 
Cm ' 



(123) 



where Eq. H22(l has been used. Recalling H119|l and making use of the steady 
intensity Eq. (|12()(l . the function tp (Cm) reads 



V'(Cm) = 



7A 



2(1+7A) 



ln7^^ 



X In 



2 + 7A 
" 2(1 + 7A) 

7|ln7^2|(A- 1) 



(I + 7A) (A + 7)+7^^^(^ 



(124) 



Equation (|123|l is the characteristic equation governing the eigenvalues A. It is a 
transcendental equation in A and hence no explicit expressions can be obtained. 
However in the class B limit 7^0, the use of an asymptotic expansion of (|123|l 
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in powers of 7 similar to that introduced in Sect. 12.2.31 allows to completely 
unfold the problem. After writing A = Aq + 7A1 + 7^A2 + • • • one obtains (SU] 



Ao = —ia , 
Ai = iaa 

ReAa = --^ i a"" - Sa'-" {A - I) + 2 (A - 1) 



1-^{A-1) 

,4 o„2 



A- 1 1 + 7^^, 



ImA2 = — (T^a 



4 

l--iA-lf 



where a — na. The RNGHI occurs at ReA2 = 0, which defines the following 
boundaries 



3 - 1) ± V - 1) 2? - 8 (A - 1) 



V 



1 + 7^2 



(125) 
(126) 



Notice the formal similarity with Eq. H94II corresponding to the UFL (for TZ ^ 1, 
V ^ 1 and then Eq. (|M|l is recovered). Now the two curves join at the critical 
point 



„ 12 

etc = 



(127) 
(128) 



The boundaries of the unstable domain depend on the reflectivity TZ through the 
function T). In the UFL T = 1 — Ti? {) [T denotes the fraction of intracavity 
power lost per roundtrip), V can be approximated as I? = 1 — 2 (T^ +T^) /3 + 
O (T'') and therefore the expressions for Ac and become 

Ac = 9+yr2(l+T)+0(T4), 

al = i2 + 8r2(i + r) + c'(r'*). 

Both quantities reduce for T = to the ones found in the UFL, namely Ac — ^ 
and al = 12. Clearly Eq. H125|) is meaningful only if the argument of the 
square root is non-negative, and this requires 2? > or, equivalently, that the 
reflectivity must be greater than 



7^n 



0.5379. 



(129) 



This is an outstanding result as it sets an absolute minimum value to TZ in 
order to observe the RNGHI (see also (SEl)- In other words, no matter how 
large the pump or how long the cavity be, the laser will not display multimode 
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Figure 6: Behavior of the critical values Ac and ac as functions of the reflectivity 



instabilities if 7^ < T^min- Above this reflectivity value, the effect of lowering TZ 
from the UFL {TZ — > 1) is monotonous: Fig. |H| shows that the critical values 
{Ac, Uc) grow by decreasing TZ and diverge by approaching 7?.niin- Accordingly, 
as shown in Fig. [7| the instability tongue H125|) shifts towards larger pumps and 
larger wavenumbers. 

To conclude, let us remark that unlike in the UFL now the value of cavity 
loss influences the second-to-first threshold ratio through TZ. In the UFL cavity 
loss affects the pump values necessary for lasing and for instability, but not their 
ratio. This fact marks a way for separating experimentally the two thresholds 
in erbium-doped fiber lasers (remember that for these lasers Wc/Won = 1 + e in 
the UFL, i.e. for small cavity losses, rendering the observation of the instability 
difficult). This is the strategy followed in the experiments IHEUHHIj which are 
treated in detail in Sect. 0] 

2.3.2 Role of distributed losses 

Another factor present in almost any experiment and not considered by the 
standard RNGH theory is the existence of distributed losses along the ampli- 
fying medium, apart from those localized outside it. The two types of losses 
should have different consequences as the role of distributed loss should be to 
decrease gain, while localized loss determines the effective reflectivity value TZ, 
whose influence has been studied above. An analysis of the combined effects of 
distributed and localized losses can be found in Ref. jlSl) which we summarize 
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Figure 7: Unstable domains for five different values of the reflectivity TZ. From 
left to right 7^ = 0.99 (which corresponds to the UFL), 0.9, 0.8, 0.7, 0.6. 

here. The proper model reads now, Eqs. lfTH|) - (pn|l . 

idr + dc)F = a{AP-xF), (130) 
drP = "/-^{FD-P), (131) 
drD = 7[l-L>-Re(F*P)] , (132) 

where we have ignored again detuning, with the boundary condition 

FiO,T)^nF{Ca,T), Cm = 27r/5. (133) 

The parameter x — ctmLm/ (|ln7?.^| + amim ), Eq. lO, verifies < x < 1, 
and am is an intensity loss coefficient per unit length ^. 

The intensity of the singlemode solution along C, and its value at the medium 
exit plane are given, respectively, by Eq. ifS^ and Eq. l(5^ with A = 0, and 
the lasing threshold is A = 1. The linear stability analysis of the singlemode 
solution follows the lines of the previous sections. As shown in [7H| the main 
effect of distributed loss is to favor the multimode instability, i.e., to reduce the 
instability to lasing threshold ratio (of course, the larger the losses the larger 
the pump needed for reaching both the lasing and instability thresholds). For 
instance, the minimum reflectivity value allowing the existence of the RNGHl, 

^We note that in Rof. 78i the parameter x was denoted by rj. Another relevant parameter 
used there, 7(j, can be written as 7d = ^ I In'R? I ^ = ^ctmim- The pump parameter A 
keeps the same definition as in 1781 . 
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Figure 8: Dependence of the minimum reflectivity necessary for the appearance 
of the RNGHI as a function of distributed losses expressed in dB. 



7?.inin7 is determined by the equation 



1 - exp {-- 



^ |ln7^2 

-X I 



X(8-9x) |ln7e2| i_ 7^2 |ln7^2 



(134) 



which can be simplified to 



1-71^ 1 



8 — 9% amim 1 — Ti-^e 



= 1. 



(135) 



This equation yields T^-min — 0.5379 for ^ = 0, Eq. p29|l . as it should. For 
nonzero distributed loss, T^-min is found to decrease with increasing ayj^L^ so 
that for ani-^m ~ 13.6, T^-min — 0.1, and for larger values of amim, T^raiu — > 
exponentially, as shown in Fig. |H1 In that figure distributed losses are expressed 
in decibels (dB) through the relation 



DL = 8.646 amir 



(136) 



We can conclude that for DL < 1 dB (which is the usual situation), T^-min 
changes very little, and only for DL > 10 dB 7?.min decreases very quickly. 
However, 10 dB is an enormous amount of distributed loss, hardly realistic for 
a usual laser. So it can safely be said that distributed losses for realistic laser 
parameters do not significantly affect the multimode emission threshold. An 
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Figure 9: Dependence on distributed losses of the minimum reflectivity neces- 
sary for the RNGHI to exist for three-level lasers. The numbers indicate the 
value of Go and the dashed line corresponds to Go — > oo. 

important question is whether this conclusion can be extended to three- and 
four-level lasers: The answer is that it cannot. For example, Eq. (|135|l is only 
valid for two-level lasers, i.e., it cannot be directly translated for three- and 
four-level lasers (this is due to the way it is derived, see the discussion in |78)'). 
When 7?.niin is calculated for three-level lasers, the result is that of Fig. |^ which 
tends to the two-level laser result only for very large G. 

In general, although the effect of distributed losses can be neglected for 
two-level lasers, as discussed above, this will not be so for three- and four-level 
lasers. The origin of this different behavior lays in the fact that distributed losses 
enter into the definition of the pump parameter, Eq. H13|l . through G, Eq. (|3()|l . 
Nevertheless, the analysis we carried out for two-level lasers leads to a simplified 
way of treating distributed loss, as long as it is not too large (roughly, smaller 
than 10 dB): Perform the linear stability analysis without taking distributed loss 
into account (which is accurate for two-level lasers), and then take distributed 
loss into account when translating the results derived for two-level lasers into 
three- and four-level lasers through G, Eq. (|30|l . We have checked that this is 
a very accurate way of deriving the stability properties for three- and four-level 
lasers for not too large distributed loss. We shall make use of this simplified 
way of treating distributed losses in Sect. I3.3l when inhomogeneous broadening 
be included 
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2.3.3 Rate equations 

Class-B lasers are usually treated with rate equations, which are derived by 
equating to zero the time derivative of the medium polarization. If this is done 
in Eqs. (|18H20I) . then the RNGHI is lost. This means that in order to capture 
the RNGH instability, coherent effects {i.e., those due to the atomic polariza- 
tion) must be preserved in the equations: The standard adiabatic elimination 
technique fails in this case. 

Nevertheless, the above circumstance does not prevent the derivation of gen- 
eralized rate equations that maintain enough information on the medium po- 
larization dynamics as for correctly describing the RNGHI. This is what we do 
in this section: We derive a reduced set of rate equations, in which the po- 
larization variable is adiabatically eliminated in a generalized sense The 
derivation of these equations stems from the observation that the RNGHI in- 
stability can be captured only if in the stability analysis terms up to order 7^ 
are kept. Therefore, the instability certainly disappears if the rate equations 
are derived by means of a standard adiabatic elimination of the polarization, 
which amounts to truncate the polarization to the term of order 7*^. But with 
a more refined adiabatic elimination, where terms up to order 7^ are kept in 
the polarization, we can derive a set of generalized rate equations that describe 
perfectly the RNGHI. 

To derive the rate equations we start from Eqs. l|^5H7n|) introduced in Sect. 
l2.1.5l for the primed variables F' , P' and D' . The most noticeable property of 
those equations is that, although they are valid outside the UFL, they admit as 
a stationary solution the singlemode solution of the UFL (|62|l . In the following 
we will consider power expansions of the dynamical variables in terms of the 
small parameter 7. To this aim it is convenient to work with variables which, 
at the leading order, are of order 1. Therefore, we introduce the new dynamical 
variables /, p and d defined as 

F' (C t) AP' iC t) 

d{Cr)^AD'{Cr). (137) 

The Maxwell-Bloch equations for /, p and d obtained from Eqs. H68|l - (|7()|l are 

{dr + dc)f = a77(C)[p-(l-*A)/] , (138) 
drP = 7"M/d-(l + *A)p] , (139) 
drd = 7 {l-t-A2-d-t-/,(C)[l- Re (rp)]} , (140) 

The new electric field / obeys the same periodic boundary condition as F' 

fiO,r) = fiCn.,r). (141) 

and in the stationary state it is not only uniform but equal to 1. In fact, Eqs. 
P38|) - (|140|) admit the stationary solution fs = 1, = 1-iA, ds = 1 + A^. For 
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a perfect resonant laser (A = 0) they reduce to 

{dr + dc)f = TviCKp-f), (142) 

drP = l-\fd-p), (143) 

drd = 7{l-d + /,(C)[l-Re(/»]} , (144) 



and the stationary solution is fs — Ps — ds = 1. To obtain the correct rate 
equations we first observe that, by setting 7 = in Eqs. (|142H144l) . we obtain: 
(i) drd = 0, hence d(C,T) = do (C), (ii) P = do (C) /, and (iii) (5, + dc)f = 
(Tri{() [do (C) — 1] /. Note that these equations do not fix the value of do (C)- 
However, since it is independent of time, it must be compatible, in particular, 
with the stationary solution d = 1, which implies do (C) = 1- Hence, at order 
in 7 we have d — 1 and f = p, from which it follows that drf = —dc^f. For 
7 7^ we make the following ansatz 

p = / + 0(7), 1 + 0(7), drf = -d^f + 0{^). (145) 

We now observe that Eq. (|143|l can be written as p = fd — jdrP, which, iterated 
twice, yields 

P = fd- ^drifd) + ^'dUfd) + 0(7') • (146) 
Taking into account Eqs. (|142|) . H144|) and (|145|) we have 

drifd) - -dac/ + ar,(p-/)-7/s(|/p-l)/ + 0(7'), (147) 
dUfd) = Oy + Oij), (148) 

which, inserted into Eq. (|146|l . yield 

p - /d + 7 H(l - d)/ + (d - j<7r^)dcf] + 7' [/s(|/|' - 1)/ + dy] + 0(7'). 

(149) 

Substituting this into Eqs. (|142|l and (|144|l we finally get the generalized com- 
plex rate equations 

idr + dc)f = (77? [/(d-l)(l- 7(777) +7(d-7fT?7)9c/] 

+7^(777 [/, (l/p - 1) / + 92^/] , (150) 

drd = 7[l-d + /s(l-|/l'd)] -^7'/sac|/|^ (151) 

which are exact up to order 7^. In the UFL we can set 77 = 1 and Is = A — 1, 
and all the parameters in the right hand sides of the rate equations become 
constant. 

As a test of this reduced model we have determined the linear stability 
of the steady lasing solution (/ l,d= 1) versus perturbations of the type 
fe(C,r) = ^a;o(C)exp(AT), with x = /, d, and ^/o(0) = (5/o(Cm), Eq. 
Solving the characteristic polynomial perturbatively up to order 7^ one recovers 
Eq. (|125|l in Sect. 12.3.11 as it must be. Notice finally that by making the limit 
7 ^ in Eqs. H150I151() (i.e., removing the terms multiplied by 7 in the field 
equation and the terms multiplied by 7^ in the inversion equation), the standard 
rate equations are recovered. 
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2.3.4 The pulse equation 

The simplest multimode solution which arises from the RNGHI consists in reg- 
ular pulses travelling along the cavity. If the homogeneous stationary solution 
is first destabilized by the iV-th sidemode, N identical pulses are present simul- 
taneously in the cavity. A semi-analytical study of the pulse solution can be 
already found in In 1989 Fu analyzed in detail the pulses and their stability 
in class-B lasers [HI]. The results of Fu were obtained in the UFL. Using the 
rate equations H15U|) and H151|) derived in the previous section we can not only 
quickly recover them, but extend them to a class-B laser outside the UFL. 

We look for a solution of Eqs. H15U|) and (|151|l which travels at velocity v 
along the cavity. To this aim we introduce the new variable r' = t — C,/v, and 
we consider the following expansions 

/ = /o(r') + 7/1 (r') + 0(7') , d = 1 + 7/s(C)rfi(T') + 0(7') , 

- = 1-7(7/3 + 0(72). (152) 

V 

We also assume that /o and /i are real and that the periodicity condition 

x{t' + NT) = x{t') , (153) 

holds for X — /o, /i, di, where T is the duration of the single pulse, and N 
is the number of pulses circulating in the cavity. Denoting differentiation with 
respect to r' by a dot, the leading order pulse equations are 

{v + P)fo-IsVfodi = 0, (154) 
di-l + /o - 0. (155) 

We can get rid of the dependence on C of /g and rj = A/ {1 + 1^) hy averaging 
both sides of Eq. H154I) on ( from to Cm. We denote this average by (...). 
Taking into account that the variation of /g along ( is ruled by Eq. H48|) . with 
X = A = 0, we find that 

(7?) = 1 , (/gry) = A - 1 . (156) 
If we define Iq = /q , the pulse equations read 

lo = \^^^ Iodi, (157) 

di = 1-/0. (158) 

The parameter (3 can be determined by means of the solvability condition of the 
next order problem, which imposes that must take one of the two values 

2(^-1) , X 
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where a± {A, TV) are the boundaries of the stabiUty domain of the homogeneous 
solution given by Eqs. H125H126ll in Sect. 12.3.11 The final form of the pulse 
equations is then 

/o - allodi, (160) 
di = 1-/0. (161) 

These equations allow to calculate with a very good approximation the pulse 
shape in the limit 7-^-0, by-passing the problem of the very long transients 
associated with the full set of Maxwell-Bloch equations, and also with the gen- 
eralized rate equations, in that limit. We have just to determine, for any pump 
A, reflectivity TZ, and cavity length = 27r/a, the correct initial conditions 
which give a periodic solution. 

To this aim we observe that Eqs. HI6UHI6II1 have a Hamiltonian structure, 
which can be better put in evidence with the change of variables x = In/p, 
p = di |31j . The equations for x and p read 

X = a^p, (162) 
p = 1-e^, (163) 

which are the equations of an anharmonic oscillator of mass subjected to 
the potential V{x) = — x. A constant of motion for this system is the energy 

Hix,p)^al^ + V{x), (164) 

from which Eqs. H162lll63'jl can be derived through Hamilton's equations. The 
anharmonic oscillator bounces between the point Xmm and Xn-mx where p — 0. 
Hence, we have 

H = V{xnun) = V{x,^,^) , (165) 
which establishes a relation between a^min and a^max- From Eq. (|164|l we obtain 

/2 

P^—^H-Vix), (166) 
a± 

which can be inserted in Eq. (|162|l to obtain the time spent by the oscillator to 
move from point xi to point X2 

1 dx 
ri2 = -7^ / , . (167) 



The period of a complete pulse is twice the time spent to go from Xmin to x 
For the A^-pulse solution it must be equal to 2tt /aN, which means 



max- 



Since Xmin and Xmax are related by Eq. (|165|l . for any given un and TZ Eq. 
p68(l allows to determine A as a function of .Tmin = hi/min (or a;max = In /max) 



40 




A 

Figure 10: Maximum intensity of the 1 -pulse solution according with Eqs. (|165|l 
and Hl()8|) for TZ = 0.7 and a — 2tt. The bifurcation is supercritical. The squares 
represent the value obtained by numerical integration of the laser equations. 

and then /min (or /max) as a function of A. Eq. (|168|l has two solutions 
corresponding to the two signs of a± . It can be easily demonstrated that the 
integral in Eq. H168|l tends to \/2'it in the limit Xminja^max 0, i.e. when 
the pulse solution reduces to the homogeneous one I — 1 [HI]. Hence, in that 
Hmit Eq. (|168|l reduces to a± — un- This means that the pulse solutions 
associated with a+ and a- are the multimode solutions emerging, respectively, 
from the upper (a+) and lower (a-) instability domain of the homogeneous 
solution. General arguments of bifurcation theory corroborated by the results 
of the numerical simulations indicate that only the former is stable. 

The plots of /max = /max(^) show that the bifurcation from the homogeneous 
solution to the pulse solution is supercritical or subcritical depending on whether 
ajv is larger or smaller than a^. In Figs. llOlandlTTlwe considered a rather small 
value of the reflectivity, TZ — 0.7. This choice of TZ corresponds to a very bad 
resonator, in which the fraction of power lost per roundtrip is (1 — TZ^) = 0.51. 
Thus, we are very far from the UFL. For this value of TZ the critical and minimum 
values of a are, respectively, Uc = 4.252 and amin — 4.212, and the critical pump 
is Ac = 13.055. In Fig. E|we consider the value aN — 27r, larger than etc, while 
in Fig. 111! we set un — 4.22, which is intermediate between amin and etc. The 
solid (dotted) lines correspond to the stable (unstable) pulse solution associated 
with a+ (q;_). When the bifurcation is subcritical, the pulse solution is stable 
also below the instability threshold up to Ac, and in this interval it coexists 



41 



2.5 




max 



2.0 - 



1.5 



1.0 



13.1 



13.2 



13.3 



A 



Figure 11: Same as Fig. ^|for a = 4.22. Here the bifurcation is subcritical. 

with the stable homogeneous solution, see Sect. 12.3.61 However, we notice that 
this happens only in the narrow range of values of ajv between amin and ac- 

In Figs. 1101 and [TTl the symbols indicate the peak intensity of the pulse ob- 
tained by the numerical integration of the dynamical equations. The agreement 
with the curve /max(^) calculated using Eqs. (|165|l and Hlt)8|l is very good. 

The stability of the iV-pulse with N > 2 was analyzed in in the UFL 
using the Floquet method. It was found that the iV-pulse solution becomes 
unstable beyond a certain value of /„iax which depends only on N. For instance 
I^^^ = 4.914 for iV = 2 and /^ax = 2.339 for = 3. For iV = 1 the stabil- 
ity analysis is more difficult and, to our knowledge, it remains an unresolved 
problem. 

2.3.5 Numerical results 

In the previous section we have illustrated a method to calculate the A^-pulse 
multimode solution developed by the laser beyond the RNGHI threshold. In this 
way we can determine at once the long-term regime of the laser, skipping all 
the transient evolution. This is an enormous advantage from the computational 
point of view, because we know that the unstable eigenvalue is of order 7^, 
which means that the time scale of the transient is of order 7"^. 

But when the pulse solution itself destabilizes, the only way to study the laser 
dynamics is the numerical integration of the dynamical equations H18l - I20|l , with 
the boundary condition H21|) ^. These are Partial Differential Equations (PDEs) 

*The pulsing regime has been studied numerically by few authors 1311 1321 1421 1431 |46I I6bl 
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that can be integrated using the finite difference method originally proposed by 
Risken and Nummedal [Hi]- However, the efficiency of this method is strongly 
limited by the the fact that it imposes the constraint of equal time and space 
step. Since the time step is proportional to 7, in the limit of very small 7 the 
computational time diverges, because the above mentioned constraint implies 
that the number of points of the spatial grid is proportional to 7^^. 

An alternative method consists in a modal expansion of the electric field 
in Fourier modes, which allows to convert the PDEs into a set of Ordinary 
Differential Equations (ODEs). The big advantage of ODEs is that they can be 
integrated using standard Runge-Kutta methods, which are easier to implement 
and generally run faster than the finite difference methods used for integrating 
PDEs. 

It is commonly believed that modal expansion is convenient only in the 
UFL, because only in that limit the electric field can be described properly by 
a limited number of modes. Outside that limit, even the stationary solution 
requires a high number modes, because it has a nontrivial spatial dependence 
along the longitudinal direction |55) . and the situation certainly worsens when 
the RNGHI sets on. 

However, it was demonstrated that in class B lasers a modal expansion of 
the electric field is fully justified for the study of the RNGHI, even outside the 
UFL, i.e. for any mirror reflectivity |8(J) . This result is a simple consequence 
of the fact that, according to Eq. (|122|l derived in Section F2. 3. II the unstable 
modes in the RNGHI depend on C as 

F,(C)exp[-A„C-V^n(C)] , (169) 

where 4'n{C) is given by Eq. I|124|l with A = A„, and An is a solution of the 
characteristic equation l|123(l . The important point is that An ~ —ia-m where 
an is the spatial frequency of the 71-th mode, and ipniO ~ in two limiting 
cases: (i) in the UFL, and (ii) in the class-B limit 7 — + 0. Then, in both limits 
at the leading order the unstable modes can be approximated by 

F,(C)exp(ianC) . (170) 

In the UFL, where Fs = \/ A — 1 can be assumed independent of the unstable 
modes immediately coincide with the empty cavity modes exp (ianC)- But this 
is true also for class-B lasers beyond the UFL, provided the scaled electric field 
/ defined by Eq. (|137|l . which is equal to 1 in the stationary state, is considered 
in place of F. Therefore, the modal expansion can be applied to the Maxwell- 
Bloch equations 1113811141)1) written in Sect. 12.3.31 for the scaled variables /, p 
and d and to the generalized rate equations that we derived from them. In both 
cases the periodic boundary condition 

/(0,t) = / (Cm , r) , Cm = 2^/5 , (171) 

and these studies have been quite superficial in most cases. 
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Figure 12: Enlargement of the steady pulse around the maximum for 7 = 0.01. 
The four outputs obtained with the four different methods described in the text 
are compared. 

allows for a modal expansion of the electric field / 

+00 

/(C,r)= e^^-VnW, (172) 

n— — 00 

and the ODEs are simply obtained by inserting this expansion in the equation 
for the electric field. For the complete Maxwell-Bloch equations the equation 
for the generic mode amplitude /„ is 

dfn Z"^" 

= -ianfn + cr y rfCr?(C) exp (-ia„C) ip- f)- (173) 

The integral can be evaluated over a grid of points Q, I = 1, . . . L. To do that 
we need to know only the values of p and d in those points. Thus, we have 
to consider a set of ODEs for the mode amplitudes /n(T) and the variables 
Pii'r) =p{Q,t) and di{T) = d{Ci,T). 

For the rate equations we proceed similarly, but with two big advantages: 
the number of ODEs is smaller because the polarization has been adiabatically 
eliminated, and the stiffness of the equations has been reduced because the fast 
time scale associated with the polarization has disappeared. 

Summarizing, we can then use three different methods to study numerically 
the laser equations: (i) the finite difference method applied to the full Maxwell- 
Bloch equations, (ii) the modal expansion technique applied to the Maxwell- 
Bloch equations, and (iii) the modal expansion technique applied to the rate 
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Figure 13: Same as Fig. E|for 7 = 0.001. Here only three curves are shown 
because the integration of the PDEs was too time consuming. 

equations. Method (iii) is expected to be more efficient than method (ii), and 
method (ii) more efficient than method (i). On the other hand, only method (i) 
is applied to the full set of Maxwell-Bloch equations, while with methods (ii) 
and (iii) some approximations are introduced. 

In Figs. 1101 and im the symbols refer to the calculations made with method 
(ii), but on that scale the results of the three methods are hardly distinguishable. 
A better check of the correctness of the three methods and of their respective 
efficiency is the comparison among the shape of the calculated pulses and that 
of the semi-analytic pulse solution found in the previous section, which is the 
result to which any dynamical simulation must tend in the limit 7 — s- 0. 

We chose to make the comparison for TZ = 0.7 and A = 30, and we consider 
the iV = 1 pulse, setting = ai — a = 27r, which implies a ~ 0.357. The 
results of the numerical simulations are shown in Fig. 1121 for 7 = 10^^ and in 
Fig. El for 7 = 10"'^. In both figures we show an enlargement of the portion of 
the pulse around the maximum to emphasize the differences among the various 
integration methods. The solid lines represent the semi-analytic solution, the 
dotted lines the integration of the PDEs, the dashed lines the integration of 
the ODEs derived from the complete model, and the dotted-dashed lines the 
integration of the ODEs derived from the rate equations. 

As expected, for the smaller value of 7 the agreement between the numerical 
simulations and the analytic result is better. As for the computation times, for 
7 — 0.01 the ratios among the three methods are 5.8:1.6:1, so their efficiency is 
still comparable. But for 7 = 0.001 the ODEs derived from the rate equations 
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are about 7 times faster than those derived from the complete model, and the 
PDEs are already too slow to converge in a reasonable time. These results 
demonstrate that to simulate the dynamical behavior of a fibre laser, for which 
7 ~ 10~^, the ODEs derived from the rate equations are the only practical tool. 

2.3.6 Subcriticality of the RNGH 

An important issue of the RNGHI that we have referred to only superficially in 
Sect. 12.3.41 is the subcritical or supercritical character of the bifurcation leading 
from singlemode to multimode emission. This issue has been considered in the 
past, as we detail below, within the UFL and may be important for the correct 
understanding of the experimental results, see Sec 01 We will here present 
some numerical results outside that limit, and therefore treat this question in 
this section. 

In their second paper of 1968 jSJi Risken and Nummedal already treated 
this question. They commented on their numerical finding that multimode 
emission persisted for a given pump even when the spatial frequency a was 
decreased below a_, see Eq. (I125|l . where the linear stability analysis predicts 
stable singlemode emission. They explicitly comment that for a < there is 
bistability between the singlemode and the multimode solutions. 

Later on, Haken and Ohno obtained an equation for the critical mode 

and found this coexistence of solutions again. In fact they derived a condition 
for supercritical or subcritical bifurcation but they did not give the dependence 
on the wavenumber a of their result because of the excessive complexity of their 
equations. Lugiato et al. j46l derived an approximated multimode solution 
(from a five-mode truncation of the Maxwell-Bloch equations) that explicitly 
showed the bistable behavior. The question was taken up again by Fu 62 who 
derived the unambiguous condition in the limit of class-B lasers that we have 
already found in Sect. 12.3.41 The bifurcation is supercritical (subcritical) when 
a > Uc (a < Qfc). Finally, Carr and Erneux j64| analytically obtained the same 
result in a slightly different limit for class-B lasers (7||/k ^ 1, which in our 
notation corresponds to small 7 and large a). 

As stated, all these studies were carried out within the UFL. The meaning of 
these results is that multimode emission could be found for parameter settings 
for which the singlemode solution is still stable. Nevertheless, in all previous 
studies, the minimum instability threshold pump Tc = 9 for class-B lasers was 
found to be a lower bound for multimode emission. Outside the UFL this ques- 
tion has not been treated before, to the best of our knowledge. The analytical 
treatment has in fact been presented in Sect. 12.3.41 where we have commented 
on this. Here we present preliminary numerical results |92| obtained with the 
the modal expansion technique applied to the rate equations (11501151(1 . In Fig. 
1141 we represent both the linear stability multimode emission threshold and the 
(numerically determined) domain of existence of stable multimode emission as a 
function of the reflectivity of the cavity for a = adTZ ^ I) = \/T2 and 7 = 0.01. 
The remarkable result is that multimode emission extends well beyond the limits 
in reflectivity marked by the linear stability analysis, although the pump nec- 
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Figure 14: The dashed line indicates the hniit of the domain of existence of 
pulsing solutions in the {TZ,A) plane for a = \/T2 and 7 = 0.01, as obtained 
by the numerical integration of the rate equations model, Eqs. (|150I151() . The 
solid line is the prediction of the linear stability analysis, Eq. (|125|l . 



essary is always larger than Ac = 9, the value theoretically predicted in Sect. 
12.3.41 This means that the domain of bistability between the singlemode and 
the multimode solutions is very large and, consequently, multimode pulsations 
could be experimentally observed for parameters for which the theory predicts 
singlemode emission. We shall come back to this later on when discussing the 
experimental observations. 

3 MULTILONGITUDINAL MODE EMISSION 
IN INHOMOGENEOUSLY BROADENED RING 
LASERS 

Up to now we have only considered ring lasers with a homogeneously broadened 
active medium. Nevertheless most lasers present some kind of inhomogeneous 
broadening in the gain line: In gas lasers the Maxwellian distribution of atomic 
velocities introduces Doppler broadening, and in solid state and fibre lasers 
differences in the sites of excited atoms introduce some spread in their transition 
frequency HH ESI El ■ 

From the point of view of multimode emission, inhomogeneous broadening 
in ring cavity lasers is of primary importance as it is well known that these 
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lasers easily emit in several longitudinal modes because of spectral hole burning 
[niE2E3lll- This mechanism for multimode emission is intuitively under- 
stood by imagining the different cavity modes, whose frequencies fall within the 
inhomogeneous linewidth, interacting independently with different subsets of 
atoms whose frequencies are quasiresonant with that of each longitudinal mode. 
If the different mode frequencies are well separated, as compared with the ho- 
mogeneous linewidth, the above is a good picture of what actually happens in 
most inhomogeneously broadened ring lasers. When the frequencies of adja- 
cent modes are moved closer together, or when their intensities grow and they 
undergo power broadening, some interaction between the modes may appear, 
mainly through saturation effects. All this has been well known for a long time 
and is not particularly exciting from the viewpoint of laser instabilities. 

Nevertheless this image may be too naive: In the previous sections we have 
analyzed in detail how Rabi-splitting induced gain leads to multimode emission 
in homogeneously broadened lasers. Then natural questions arise immediately: 
How does this mechanism affect multimode emission in inhomogeneously broad- 
ened lasers? Or in other terms, how does inhomogeneous broadening affect the 
RNGHI? Are there two different mechanisms, namely Rabi-splitting induced 
gain and spectral hole burning, for multimode emission in inhomogeneously 
broadened lasers? 

The rest of this section is organized as follows. In Sect. 13.11 we present the 
Maxwell-Bloch and rate equations models. Then we analyze the multimode 
emission threshold, first in the UFL (Sect. I3.2|l and then outside this limit 
(Sect. 13 .311 . In Sect. 13.21 we pay special attention to the comparison between 
Maxwell-Bloch equations and rate equations. 

3.1 Modelling 

As in the previous sections, we consider an incoherently pumped two-level ac- 
tive medium of length L^, contained in a ring cavity of length Lc, interacting 
with a unidirectional plane wave laser field. However, now the field interacts 
with an inhomogeneously broadened active medium. We model the medium by 
considering a Lorentzian distribution of atomic frequencies. We further assume, 
for the sake of simplicity, that the cavity is resonant with the center of the 
atomic transition frequency distribution. 

3.1.1 Maxwell Bloch equations 

The Maxwell-Bloch equations describing such a laser can be easily written by 
generalizing Eqs. H18II20|I . The generalization consists in replacing the polar- 
ization appearing in the field equation Eq. H18|l by the weighted sum of the 
polarizations corresponding to the different frequencies lo. The model equations 
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read [El EH E3| 



{dr + dOFiCr) 



aA 



dL0C{L0)P , 



(174) 



a,P(u;,C,T) 
a,D(u;,C,T) 




-{l + iuj)P + FD] , 
D-Re{F*P)], 



(175) 
(176) 



supplemented by the boundary condition 



F{Q,T)=nF{Cra,T) , 



C. 



27r/a. 



(177) 



Notice that distributed losses have been neglected (x — 0) and resonance has 
been assumed between the cavity and the center of the atomic frequency distri- 
bution ((5 = 0). 

In Eqs. (|174|l - H17t)|l F{C,, r) is the normalized slowly varying envelope of 
the laser field, and P(w, C,t) and D{u},(^,t) are the normalized slowly varying 
envelopes of the medium polarization and the population inversion, respectively, 
for molecules detuned by oj with respect to the cavity resonance. 

Inhomogeneous broadening is accounted for by the function /3(w). It repre- 
sents the spectral distribution of atomic resonances, which, in order to deal with 
analytical expressions, is taken to be a Lorentzian distribution of width (HWHM) 
u 



With our notation the unsealed total gain linewidth (HWHM) of the medium 

is 7_L (1 + u). 

3.1.2 The uniform field limit 

As demonstrated in Sect. 12.21 when the effective amplitude reflectivity TZ is 
close to unity, one can apply the uniform field limit that consists in replacing 
the field Eq. ((TTljl by 



complemented with the new boundary condition _F(0,t) = P(<^,ii,t). 

3.1.3 Standard rate equations (uniform field limit) 

Class-B lasers (those for which 7j_ 7||, k, i.e., 7, c <C 1), use to be described 
with rate equations, which are obtained after the adiabatic elimination of the 
medium polarization in the Maxwell-Bloch equations. These rate equations are 
generally assumed to be appropriate for describing multimode emission due to 
spectral hole burning ^Jj. As we shall later compare the predictions of the full 
set of Maxwell-Bloch equations with those of the simpler rate equations, we 





(179) 
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give here this simpler model, whose detailed derivation can be found in |11[I94| . 
After expanding the field as 



+N 



F(C,t)= /^e^^" exp [za„ (C - r)] 



(180) 



-N 



with a„ = na, expanding the medium polarization in a similar way and adiabat- 
ically eliminating the polarization, one is left with the standard rate equations 



din 
dT 



= 2cr/„ 



A 



duj 



auj)D 



dD 
97 



7 — 7Z3 

The phases can be calculated through 

-uA 



'-Si 



1 + (w - a„)2 



1 



+ (ti; - UnY 



dT 



C{lu)D{lu - an) 

UiO— 7 rij— 

-00 1 + (w - anY 



(181) 
(182) 

(183) 



In these equations, the integer index n denotes the n-th longitudinal mode with 
respect to the central resonant mode (n = 0). As we show below, the rate 
equations description is not appropriate for too long cavities, where long may 
actually be quite short. 



3.2 Multimode emission in the uniform field limit 

The stability properties of the uniform field limit were treated first by Mandel 
inni IMl (see also Q^D; who showed that in the limit u — > 00, the multimode 
emission threshold is just above the lasing threshold. The general case has been 
treated at depth in [^, and the reader is referred to that publication for full 
details. 

Eqs. (| 17511 761 179() have two stationary solutions: The laser-off solution 
(F = 0, P = 0, D = 1), and the monomode solution, for which the intensity of 
the lasing mode, /q = F^, can be written as ^ = \/l + /q (Vl + Iq + u). The 
threshold for lasing emission is found by taking /o = in this expression, and 
thus it occurs at a pump A = Aq = \ + u. Then, by defining the normalized 
pump parameter r = A/Aq, the lasing solution reads 



RiR + u) 
1 + u 



(184) 



R = ^l + Io, (185) 



or 

Io = r 



(1 + It) - 1 + - y^u^ + 4r (1 + u)) , (186) 

which for u = gives /q = r — 1. At threshold, R = r = 1. 
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3.2.1 Two estimates of the multimode emission threshold 

Let us assume momentarily the following naive approach: Assume that the 
threshold for amplification of a detuned mode is not affected by the already 
existing resonant lasing mode. This assumption will be reasonably accurate 
whenever (i) the inhomogeneous width u is large, because in that case the 
intensity Iq of the resonant mode at threshold for multimode emission will be 
small, and (ii) when the frequency of the detuned mode is (sufficiently) larger 
than the normalized homogeneous width 7"^ (which is a measure of the width 
of the spectral hole), i.e., for large a. 

In [5T] we showed that under these assumptions the emission threshold for 
a sidemode of spatial frequency a can be estimated to be given by 

\ja + 1 + u J \l + u/ 

where the limit corresponds to class-B lasers (7^1). Thus multimode emission 
occurs, in this simple approach, for r > rthr- The term 7a/(l + u) represents 
the sidemode frequency offset divided by the sum of the cavity linewidth and 
the total gain linewidth. One can think of Eq. (|187|) as the multimode emis- 
sion threshold when spectral hole burning is the only relevant mechanism for 
multimode emission and there is not any interaction between modes. 

A second, more accurate, estimate for the multimode emission threshold for 
class-B lasers is obtained from the rate equations model, Eqs. (|181l 118^ . The 
derivation can be found in Chapter 4 of Khanin's book ^T] (see also (Slj)- It 
reads 

^ p{R)+-f^a^, (188) 
piR) = R'^{l-u)+R{l-u){u + 2) + {u + l)'^ , (189) 

when the cross relaxation coefficient 7cr is set equal to zero in Khanin's equa- 
tions JT]. By using Eq. H184|) . the above equation can be rewritten as 

, , 2R{R+1) 

rrate = rthr.s(a) + — ^ "2^, (190) 

(l + u) 

with rthr,B (a) given by Eq. (|187|l and R given by Eq. (|185|l . the monomode 
solution being unstable for r > rratc- We see that riatc > ?'thr,B(Q;), i.e., the 
instability threshold given by Eq. (|190|) is always larger than the one predicted 
by Eq. (|187|l . which was admittedly valid for a 3> 1. Both boundaries tend 
to coincide in the limits of validity of Eq. (|187|) . i.e., u ^ 1. This is perfectly 
reasonable as Eq. f ll9U|l contains one more ingredient with respect to Eq. H187fl . 
namely the saturation of the sideband gain due to the resonant mode, which 
explains why it gives a larger value for the instability threshold. 

Notice that Eq. (|190|l predicts that there will be no multimode emission un- 
less M > 1, i.e., it fails in predicting multimode emission due to Rabi-splitting 
induced gain (RNGHI). This is entirely normal as the rate equations do not 
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contain information about the atomic coherence, i.e., about the medium polar- 
ization dynamics, which is the responsible for the RNGHI as we have seen in 
the previous section. On the other hand, the limitation u > 1 does not appear 
in the naive approach of Eq. (|187|l . Then the physical picture is the following: 
For u < 1, even if the threshold for emission for more than one mode is crossed, 
the nonlinear competition between the different cavity modes enforces that only 
one of them lases. This is an example of "the winner takes all" competition and 
is the common picture of multimodc dynamics — which we know to be false 
because of the RNGHI. 

The existence of a minimum instability threshold in the limit a ^ pre- 
dicted by the rate equations, Eq. 1)1901) . is quite unphysical, because it would 
mean that the closer are the sidemodcs to the central mode, the lower is its 
instability threshold. This erroneous result is a consequence of the approxima- 
tions introduced in the derivation of the rate equations, and it disappears when 
the stability analysis of the homogenous solution is performed using the full 
set of Maxwell-Bloch equations. As will be shown below, the complete model 
predicts that, for a given pump value, there is always a minimum value of the 
frequency spacing a below which there is no instability because the gain of the 
sidemodes is saturated by the strong central mode. 

3.2.2 Multimode emission threshold 

The multimode emission threshold is rigorously obtained by performing a lin- 
ear stability analysis of the single-mode solution Eq. H184|) . The analytical 
expressions obtained are cumbersome and it is not trivial to extract analytical 
information. We refer the interested reader to |81j where several limits of inter- 
est are treated explicitly. Here we shall make a resume of the results concerning 
class-B lasers. Two limits of interest can be treated explicitly for class-B lasers: 
The short and long cavity limits, i.e., the limits a 3> 1 and a = 0{1) or smaller. 

For short cavities, a 1, the instability threshold turns out to be given by 
Eq. (|190|l . This amounts to saying that for short cavities the rate equations 
approach is valid for deriving the multimode emission threshold. 

For long cavities, a = 0(1), an approximate expression for the multimode 
emission threshold can be derived. It is given by 

r = 0, (191) 
V = p (R) - SiR^ ~ I) [{R + u + If - Ru] 
+R{R^ - 1){R + lf{R + u){u + 2), 

with p {R) given by Eq. (|189|l . The monomode solution is unstable for V < 0. 
Before comparing this multimode emission threshold with that predicted by the 
rate equations, we show the influence of the inhomogeneous broadening on the 
multimode instability threshold. In Fig. 1 151 we represent the multimode emission 
threshold predicted by Eq. (|191|l for several values of u in the (r, a) plane. 
Notice the enormous quantitative effect that the inhomogeneous broadening has 
on the instability threshold. If fact, it can be shown analytically FT" that the 
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Figure 15: Multimode emission threshold in the (r, a) plane for the values of 
the inhomogeneous to homeogeneous linewidth ratio, u, marked in the figure. 
The thick line corresponds to the homogeneous limit u — Q. 

critical point (rc, ac), that for which r is minimum at the instability threshold, 
varies as 



for u » 1. 

In Fig. Iinithe thresholds predicted by Eq. lfTH5|) (dashed line) and Eq. ((TnT|l 
(full line) are represented in the (r, a) plane for u = 4. The exact instability 
threshold (diamonds) as obtained by numerically solving Eq. (22) in |SJ, is also 
shown. Several conclusions can be extracted from this figure. 

(i) There is a single instability threshold, i.e., there are not two mecha- 
nisms for multimode emission (spectral hole burning and RNGHI) but a single 
mechanism; 

(ii) The asymptotic expressions Eqs. H188() and (|191|l compare perfectly well 
with the exact result in their respective domains of validity: The threshold 



9 - 36m + 0{u^) , 
12 - 33.75m + 0(^2), 
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for M ^ 1, and as 
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Figure 16: Multimode emission threshold for u = 4. The sohd (dashed) hnes 
have been obtained from Eqs. (|191|l and H188|) . respectively. The diamonds 
correspond to the exact stability analysis of the Maxwell-Bloch equations, Eq. 
(22) in Ref. [HI], and have been calculated for 7 — 10~^ and cr = 0.1. 

predicted by the rate equations becomes invalid for small a (long cavity) while 
the approximate expression Eq. p91|l becomes invalid for large a (short cavity); 

(iii) The unphysical instability predicted for small a by the rate equations 
result is removed, and there is a minimum value of a below which there is no 
instability; and 

(iv) there is a domain of intermediate cavity lengths where the two analytic 
expressions connect. 

The last item allows to determine which is the value of the cavity length 
beyond which the rate equations approach is no more valid. By analyzing Eqs. 
H188|l and (|191(l one obtains that the distance between the two functions is min- 
imum for a ~ 2j~^/^ when u > 2. This leads to a "coherence length" estimate 



for lasers with a significant inhomogeneous broadening (u > 2). 

Then, for cavities larger than icoh, atomic coherence effects {i.e., Rabi- 
splitting induced gain) are important for the determination of the multimode 
emission threshold while they are irrelevant for shorter cavities. For example, 
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for CO2 lasers Lcoh ~ 100 m, for 632.8 nm HeNe lasers Lcoh ~ 20 m, for Er^+- 
doped fibre lasers icoh ~ 5 cm, and for Nd-glass lasers Lcoh ~ 3 mm. Thus, 
the necessity of using the full Maxwell-Bloch description of the laser in order 
to describe multimode emission depends strongly on the particular laser system 
under consideration: It is necessary for Nd-glass and Er^+-doped fibre lasers 
and unnecessary for CO2 and HeNe lasers. 

Let us remark that the expression for Lcoh gives a good estimate for the 
critical length for values of m > 2; for smaller values of u, Lcoh decrease, tending 
to for u ^ 1, so that the necessity of considering coherent effects is still more 
important in this limit. 

The above shows that Rabi-induced sidemode gain (RNGHI) is far from be- 
ing a negligible mechanism for multimode emission in inhomogeneously broad- 
ened ring lasers. Quite to the contrary, it is an essential ingredient for those 
lasers with cavity lengths larger than Lcoh or even shorter if the ratio of inho- 
mogeneous to homogeneous linewidth is small. 

3.2.3 Dynamics beyond the multimode emission threshold 

The numerical integration of the inhomogeneously broadened model is far from 
being a simple task, especially for class-B lasers. In fact, this task has been 
carried out very recently for the first time in [M]. For the sake of illustration, 
let us comment that some of the calculations we comment below implied a 
running time of more than one month on a R12000 Silicon Graphics processor. 
This illustrates the enormous difficulties of making numerical simulations for 
class-B laser parameters, due to the stiffness of the equations in this limit. 

There were previous attempts, during the eighties, by Brunner et al. 
EE] of numerically studying multimode emission in inhomogeneously broadened 
lasers. But the model they used was that of a standing-wave cavity and the 
third-order Lamb theory approximation was assumed. As this approximation 
does not contain atomic coherence effects, their results did not give information 
about the RNGHI mechanism in standing-wave lasers. 

In IHH, the Maxwell-Bloch Eqs. (|175I176I179|) and the rate equations model 
Eqs. (|181I182|I were numerically integrated for fixed pump (r = 5), inhomo- 
geneous broadening (u = 2) and spatial frequency (equal to the homogeneous 
linewidth), and three different values of 7 (10~^, 10~^, and 10^^). As we change 
the value of 7 over three orders of magnitude, we change the Rabi frequency of 
the field. In doing so, the behavior makes the transition from that predicted 
by the rate-equations model (for smaller 7) to that in which Rabi-induced 
gain effects are expected to be important (for larger 7). The main finding in 
|M] consists in the confirmation of the conclusions extracted from the stability 
properties: Rate equations describe properly the multimode dynamics only for 
short cavities, that is, under conditions for which the free spectral range is large 
as compared to the Rabi frequency of the intracavity field. In this limit rate 
equations and Maxwell-Bloch equations provide the same results for the mode 
intensities at steady state, and no phase locking is found even in the Maxwell- 
Bloch equation, where all relevant information about the phase dynamics is 
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Figure 17: Multimode emission threshold in the (r, a) plane for m = 2 and the 
values of the reflectivity marked in the figure as obtained from Eq. (36) in |93| . 

preserved (remember that mode-locking cannot be described by the rate equa- 
tions model as modal phases do not enter in the dynamics of mode intensities 
and population inversion, see Eqs. H181l - I183'jl 'l. Then, for short cavities the use 
of rate equations is perfectly legitimate. 

This is no longer true for longer cavities, where the inhomogeneously broad- 
ened laser behaves like a homogeneously broadened one, with a multimode dy- 
namics which presents all the features of the classical RNGHI: The number of 
excited modes is larger than expected from linear stability considerations and 
phases spontaneously lock, giving rise to regular pulses in the total intensity. 
We refer the reader to (Ml for full details. 

3.3 Multimode emission outside the uniform field limit 

Outside the uniform field limit, i.e., for arbitrary values of the amplitude re- 
flectivity TZ the analytical expressions for the instability threshold are much 
more involved than in the uniform field limit previously analyzed. For class-B 
lasers the problem has been treated in detail in jHSJ, and the interested reader 
is referred to that paper for full details. Here we shall stress some important 
points. 

As it happens in the homogeneously broadened model, a decrease in the 
value of TZ leads to an increase of the multimode emission threshold. This can 
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Figure 18: Minimum value of the reflectivity for the multimode instabiUty to 
exist, 7?,niin, as a function of the inhomogeneous to homogeneous Hnewidth ratio. 

be seen in Fig. 1 171 where the multimode instability threshold is shown, as usual, 
in the (r, a) plane for u = 2 and several values of the amplitude reflectivity 
TZ. For Ti. = 0.99, the instability threshold is very close to that predicted 
by the homogeneously broadened model (Eq. I^^J), and decreasing values of 
TZ increase the instability threshold as expected. What is unexpected is that 
multimode emission exists for values of TZ as low as 0.1. This is unexpected 
because we showed analytically in Sect. 12. 3. II that in homogeneously broadened 
lasers, the decrease of TZ leads to the disappearance of the multimode instability 
when TZ < 7?.min — 0.538. 

Fortunately, T^-min can be analytically determined even outside the uniform 
field Hmit (see [HSl)- It is given by V {TZ, u) = with 



which reduces to Eq. (|126|l for u = ^. 

In Fig. ^] TZi-ain is shown as a function of u . Notice that the limitation 
on the value of TZ for the existence of multimode emission disappears at m = 
1 (although in practical terms it disappears for u > 0.9 as TZmin decreases 
exponentially as u — > 1). Thus the multimode instability does not have any 
limitation in cavity losses when inhomogeneous broadening is about the same 

^Notice that Eq. 11971 corrects a typographic error in Eq. (46) of inU. 
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as the inhomogeneous width or larger. This marks an important difference with 
the case of homogeneous broadening. The fact that this quahtative change 
appears at u = 1 can be related to the fact that the rate equations model 
predicts the multimode instability only for u > 1. This result suggests that 
for u > 1 spectral hole burning makes multimode emission possible when the 
Rabi-induced gain mechanism could not do it by itself (i.e., there is no more 
a "^min, which is a characteristic of the RNGHI), while for u < 1, T^-min exists 
(although its values decreases for increasing u). This suggests that spectral hole 
burning plays a less important role in this limit. 

This result is particularly relevant from the experimental point of view: An 
increase of the cavity losses eventually leads to the impossibility of multimode 
emission when u < 0.9 {i.e., for homogeneously broadened lasers or lasers with 
small inhomogeneous broadening), but this does not happen for lasers with 
enough inhomogeneous broadening at u > 0.9. 

3.3.1 Effect of distributed losses in three level and four— level lasers 

In Sect. 12.3.21 where homogeneous broadening was considered, it was shown 
that distributed losses have quantitative importance for the determination of 
the multimode emission threshold for three- and four-level lasers, although it 
was not really important for two-level lasers as far as distributed losses have 
reasonable values, roughly below lOdB. The reason lies in the way how losses 
enter in the transformation of the pump parameter for applying the results of 
the two-level theory to three- and four-level lasers. 

For inhomogeneously broadened lasers the same argument holds, and dis- 
tributed losses must be taken into account when calculating the multimode 
instability threshold for three- and four-level lasers. Performing the linear sta- 
bility analysis with inhomogeneous broadening and distributed losses outside 
the uniform field limit is a very heavy task that does not provide any analytical 
insight. Rather than going to that trouble, it suffices at least in the case of not 
too large distributed losses to use the results of the linear stability analysis out- 
side the uniform field limit without distributed losses, and introduce them when 
calculating the transformation of the pump parameter, as already explained in 
Sect. 12.3.21 This is the way we calculated the instability threshold presented 
in [Mj and reproduced in the following section (cf. Fig. 1211 below*) . For more 
details, see Ref. 24 in P|. 

4 ON EXPERIMENTAL STUDIES OF THE 
MULTIMODE INSTABILITY IN EDFL's 

In this section we describe experimental observations that contribute to the dis- 
cussion of multimode laser instabilities. We already discussed in Sect. Qlthat 
experimental research on the RNGHI has been limited to the dye laser 44, 45j 
and that the observed multimode emission was explained as being caused by 
the band structure of the lasing levels |47[ 1481 0^1 1 not a manifestation of the 
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RNGHI. Prior to Refs. [^[72] no other experimental work on the subject has 
appeared. The reason is probably that the constraints imposed on the cavity 
length by Eqs. (|93|l and H1UI)|I have discouraged experimentalists: According 
to Eq. H1UU|) . cavity lengths need to be enormously long (in some cases sev- 
eral orders of magnitude longer than standard values), while at the same time 
there is the requirement of a very large instability to lasing threshold ratio (the 
"factor-of-nine" ) By now it is known, however, that the "factor-of-nine" is 
irrelevant if the active medium is better modelled as a three-level laser medium. 
Moreover, for such lasers Eq. 1)1 UU|) does not constitute a hard limit, as for cavity 
lengths well below this value instabilities occur at accessible pump values. This 
makes it much more likely that the RNGHI becomes observable in some lasers. 

By this reasoning, Er-doped fiber lasers are rendered the most promising 
candidate for a clear observation of the RNGHI: these are three-level lasers, 
and their cavities can be made very long. We will therefore concentrate on this 
laser type. 

Many laser applications require monochromatic operation, and therefore sev- 
eral researchers attempted to operate Er-doped fiber laser in single mode op- 
eration. However, after many failures and a few partial successes, this is now 
considered as notoriously difficult. Where it was attempted, researchers chose 
one of these strategies: (i) increase the cavity's effective free spectral range 
beyond the gain bandwidth (see, e.g., [HZI), or (ii) reduce the gain bandwidth 
below the free spectral range by insertion of filters with very narrow bandwidth 
(see, e.g. HH). None of these approaches led to a full success. For example, 
in a publication subsequent to [HHl it is reported that single mode operation 
became impossible to maintain above a certain pump level |99j . The reason for 
this behavior was not further investigated. 

The pragmatic conclusion is that EDFL's obviously have a natural tendency 
to operate on several modes simultaneously. The small free spectral range caused 
by the necessarily long resonator contributes to the difficulty, but by itself cannot 
explain the underlying reason. 

4.1 A first dedicated approach 

Probably the first experimental study of fundamental causes for the instability 
of an EDFL was presented in '7T]. The authors made a distinction between 
relaxation oscillation and self mode locking, and investigated the latter. They 
used a resonator consisting of a WDM coupler to bring in the pump light, and an 
output coupler to steer fully 90% of the power out. A polarization-insensitive 
optical isolator enforced unidirectional operation. An Er-doped active fiber 
provided gain; alternatively either a 15 m length at low dopant concentration 
(300 ppm), or a 0.8 m length at high concentration («5000 ppm) was used. Al- 
lowing for some extra length of component pigtails, this brought the resonator's 
free spectral range to about 10 MHz or 70 MHz, respectively. 

The laser output was directed through another isolator and then monitored 
either by a fast photodiode hooked up to a fast oscilloscope or an RF spectrum 
analyzer, or by a background-free autocorrelator. No stable cw emission was ob- 
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served. For the longer cavity it was reported that self-modelocking produced a 
train of pulses with a repetition rate given by the resonator's free spectral range. 
An occasional presence of satellite pulses at intermediate times was mentioned. 
The pulses were reported to have a temporal width (FWHM) of a few ns. In- 
spection with the autocorrelator was performed to check for substructure. The 
autocorrelator could monitor a time window of 50 ps, and within this window 
no substructure was found. For the shorter cavity, a modulation of the output 
power at the free spectral range was also observed. However, it did not take 
the form of pulses but rather had a nearly sinusoidal shape. This is plausible 
because for a larger free spectral range, fewer modes will fall into the bandwidth 
of the gain. 

When the pump power was varied near lasing threshold, in the shorter cavity 
the instability was seen whenever there was lasing, while in the longer cavity a 
small interval of single-mode operation seemed to exist just above lasing thresh- 
old. 

4.2 The follow-up 

In a subsequent study 77, basically the same group of authors replaced the 
polarization-insensitive isolator with a polarizing isolator. They also added po- 
larization controllers to the cavity. Two output couplers branching out 95% 
or 50% were used alternatively. The active fiber was again of the low Er con- 
centration type, and was 13 m long. The cavity free spectral range was thus 
near 10 MHz. Again, self mode locking at the cavity free spectral range was 
observed; pulses had a duration of 1.74. . . 3.00 ns. The RF spectrum contained 
« 250 beat notes, indicative of a similar number of oscillating modes. 

In this experiment there was a combination of fiber birefringence plus po- 
larizing elements in the cavity. This raises the issue whether mode locking due 
to Nonlinear Polarization Rotation (NPR) might have occurred. The authors 
argue that NPR can be ruled out for the following reasons: (i) In the previous 
setup, there definitely was no picosecond structure, and when the polarizing 
isolator was introduced, the pulse shape was not modified. This suggests — 
somewhat indirectly — that the polarizing action is not responsible. The au- 
thors further noted that (ii) at times intermediate to the pulses there was a 
constant background of random signal, possibly satellite pulses of some kind, 
and that NPR would likely suppress such structure. However, random groups 
of pulses are routinely seen in NPR lasers. Finally, (iii) power levels in the fiber 
were deemed insufficient for NPR, in particular since a wide core fiber was used. 
NPR requires remarkably little power, however, and at the powers stated, NPR 
cannot be ruled out entirely. However, the authors kindly inform us that the 
instability also existed just above threshold, and in that case NPR is indeed 
highly unlikely to occur. 

There remains an unresolved discrepancy about the measured spectral shape 
of the pulses which is reported to be Gaussian and, for 2 ns pulses, must have 
been 0,001 nm wide. The optical spectrum analyzer reportedly used for this 
measurement, however, has a specified spectral resolution of 0,05 nm. 
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One of the most prominent features of RNGHI, namely its threshold-like 
onset, was not addressed in the experiment of (72,. The publication provides a 
comment that the experiments took place far above lasing thresholds, but that 
instability persisted down to at least a less-than-tenfold threshold power. As 
the authors kindly inform us, the instability was seen immediately above lasing 
threshold. From this information one must conclude that the matter deserves 
more clarification before anything is definitely proven. 

4.3 A systematic assessments of thresholds 

In subsequent work, an Er-doped fiber ring laser specifically designed for obser- 
vation of instabilities was set up by some of the present authors. As a starting 
point a standard configuration was chosen, see Fig. 1191 The cavity contained 
8.2 m of active fiber (585 ppm Er dopant level) in a 22 m long ring (the remain- 
der consisted of standard single mode fiber). Pump light came from a 100 mW, 
980 nm laser diode; it was launched into the ring by a WDM coupler. Light 
coupled out by this WDM coupler, and also from an additional 95/5 coupler, 
was used to monitor the system simultaneously by a fast photodiode and an op- 
tical spectrum analyzer. An optical isolator ensured unidirectional operation, 
and a moderately narrow bandpass filter gently restricted the bandwidth avail- 
able for lasing to about 1 nm. As a unique feature, variable loss was inserted 
into the cavity. To this end either an amplitude modulator or a tight fiber 
coil with well-defined radius was employed. High variable loss served to bring 
out the onset of instability more clearly. Great effort was made to character- 
ize the exact amount of loss for each setting: all components (localized losses) 
including splices etc. were tested individually, and the distributed loss in the 
Er fiber was determined. Finally, as a cross-check the total loss was calculated 
from the operational laser's output power vs. pump power relation. For more 
detail see [HHl- Just above lasing threshold this laser indeed operated in a 
single longitudinal mode. As the pump power was increased, modulations of 
the output power appeared at frequencies which were integer multiples of the 
cavity round trip frequency (9 MHz) without exception. Obviously these were 
beat notes between different longitudinal modes. 

These beat notes, however, were neither steady in amplitude, nor in fre- 
quency: The frequencies involved hopped rapidly and apparently at random all 
across the range up to several GHz, and during much of the time in between 
there was no beat at all. After the insertion of the bandpass filter already men- 
tioned above, the frequency hopping range was limited to hundreds of MHz and 
thus more manageable. Still, the beats remained unsteady. A typical beat note 
episode lasted on the order of tens of milliseconds to a few seconds. The tempo- 
ral profile of the beat note was almost always very nearly sinusoidal, indicative 
of a beat between only two resonator modes. In the presence of the high loss 
intentionally introduced here, the laser could not be pumped very far above 
threshold, so that it may be not too surprising that there was just dual-mode, 
but no multimode operation. 

We wish to point out that for this experiment nonlinear polarization rotation 
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Figure 19: Experimental setup. The modulator serves to introduce well-defined 
loss; it can take the form of either an electro-optic modulator or a tightly wound 
fiber coil. For further detail see text. 

(NPR) |1UI)| can be safely ruled out. Polarization-dependent losses, on which 
NPR hinges, were carefully avoided. 

Given a random phenomenon, statistical means were employed to character- 
ize it. Time series of the instantaneous power were recorded which — in view of 
the typical timescales — were chosen to be several seconds long. Unfortunately, 
adequate sampling to correctly assess all up to the highest frequencies would 
have required data rates of GB/s and file sizes of about 10 GB per shot; that 
is just not feasible. To keep data files at manageable size, undersampling at 
20 000 samples per second was chosen. While this way information about the 
actual beat frequency is lost, all episodes of mode beating longer than 100 /iS 
still can be detected from the time series. Shorter episodes seemed not to occur 
anyway. Occasional occurrences of relaxation oscillation were easily identified 
in the file by their very different amplitude, and were discarded. The fraction 
of time during which valid mode beats were detected served as the measure of 
instability M . 

The data-taking procedure consisted in setting a particular loss value, then 
incrementing the pump power in small steps, and determining M at each step. 
(At the same time, as described above, the total power was recorded to help 
assess the loss value). Next, the loss was incremented, and the procedure was 
repeated until the accessible range of loss values was exhausted. After an eval- 
uation of these extensive data, the following conclusions are reached: The laser 
power data in Fig. |2D1 show the universally known threshold behavior: the 
power is close to zero below threshold, but not exactly so due to fluorescence. 
Above threshold the power makes a good fit to a straight line by which the slope 
efficiency is defined. The measure of instability, on the other hand, strictly re- 
mains at zero up to some point above the laser threshold. Beyond that point, 
M sharply sets on to nonzero values, and continues to rise as the pump power 
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Figure 20: Laser power (filled circles) and measure of instability (open squares) 
as a function of pump power. 



is raised further. This is clear evidence that the instability has a well-defined, 
sharp onset, a fact which had not previously been demonstrated experimentally. 

In the next step both the lasing threshold and the instability onset are 
taken from data as in Fig. 1201 Fig. |^ reveals that the ratio of both values 
-Pinst / (where Pinst and Pq denote the pump power at the instability and lasing 
thresholds, respectively) does indeed scale with loss. At low loss, Pinst/-Fb ~ 1- 
Under typical operating conditions of Er fiber lasers, losses would be even lower, 
and experimenters would be unable to tell apart both onsets. For larger loss, 
however, Pinst/Po increases up to about 1.5. The intentionally high loss of this 
experiment pays off nicely here: Both onsets are clearly distinguishable, and 
the interval of single mode operation in between is clearly identifiable. This 
constitutes a considerable progress over previous work. 

However, we must emphasize an important fact regarding the range of loss 
values used in that experiment. We showed in Sect. 12.3.11 that in a homoge- 
neously broadened gain medium T^-min describes a maximum value of loss beyond 
which there exists no instability at all. 7?.min as defined in Eq. 11291 corresponds 
to 4.5 dB of localized loss in the experiment ^89,, but Fig. I21l clearlv shows that 
instability persists at much higher loss. 

To resolve this discrepancy, we reconsider the structure of the Er gain line. 
In Sect. 13.31 we showed that in the case of an inhomogeneous contribution 
to the line as expressed by u 0, T^min goes to zero (see the discussion of 
Eq. 11971 and Fig. I18() : hence the corresponding maximum loss diverges. It is 
difficult to make precise statements about the value of u for the fiber used in the 
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Figure 21: Influence of localized cavity loss. Shown are (a) the laser threshold, 
(b) the instability onset, and (c) the ratio of both. The dashed line in (a) is 
a fit with theory (see [HHl for more detail); while the line in (c) only serves to 
guide the eye. 



experiment in but a choice of u = 2 is reasonable. For this value, Fig. 
shows the theoretically expected threshold ratio as a function of localized cavity 
loss (the calculation takes distributed losses into account, see ^22, for details). 
Evidently, the disagreement with experimental data is reduced dramatically. 
While quantitatively r is systematically predicted too high, in particular at the 
highest loss values, at least the discrepancy about the existence of an instability 
onset is resolved. 

Finally we need to address a caveat about the interpretation of the exper- 
imentally observed instability onset. As was discussed in Sees. 12 .3. 412 .31)1 the 
possibility exists that the instability threshold is either a supercritical or a sub- 
critical bifurcation. Based on experimental data alone, a decision between these 
possibilities cannot be made. Consider a subcritical bifurcation, which would 
in all likelihood imply a range of bistability between the single-mode and the 
multimode solution (at least this is the case in the homogeneously broadened 
case, see Fig. I18f) . The observed intermittent behavior would make some sense 
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Figure 22: Ratio of instability onset to lasing threshold, as predicted by a 
three-level laser model with inhomogeneous broadening at m = 2. Different 
modes become unstable at different points; the labels refer them to the central 
lasing mode. For the instability onset, always the lowest line counts. 

in that case: The experimentally determined instability onset would then be the 
lower limit point of the unstable branch, and the bifurcation proper was never 
reached due to limited available pump power. In fact, not even the point of 
"Maxwell's construction" (where both branches are occupied for equal amounts 
of time on average) was reached. Even if unlimited pump power had been avail- 
able, it is not at all clear whether a final conclusion about the nature of the 
bifurcation could have been reached, because a sizeable increase of the pump 
power eventually brings on other processes like Brillouin scattering, thermal 
effects, etc., which further complicate the issue. We must therefore leave this 
question open for now. 

Let as finally remark that, in contrast to j72| . multimode operation and 
pulsing was never observed in ,;8^: there was only dual-mode operation and 
sinusoidal modulation. Also, there were only intermittent, not steady, beat 
notes between modes. Surely, this must have to do with the fact that the laser 
was intentionally made lossy so that it could never be pumped very much above 
its first threshold. Whether there are still other factors involved (codopants of 
the fiber, etc.) must remain unresolved at this point. 
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5 CONCLUSION AND OUTLOOK 



We have introduced the different models required for the study of the RNGHI. 
In particular, we have treated the applicability of the two-level laser model 
to three- and four-level lasers and, importantly, the rigorous derivation of the 
uniform field limit. We have then revised the basics of the RNGHI, and we have 
reviewed our continued work on the subject over recent years. This research 
was motivated by the suggestion by Lugiato and coworkers in 1997 that the 
pulsations exhibited by a unidirectional EDFL could be a manifestation of this 
elusive phenomenon. 

At long last, the threshold- like onset of the multimode instability in an 
EDFL was demonstrated in 89 . Data show a qualitative and, with appropriate 
corrections, even semi-quantitative agreement to theoretical expectations. Nev- 
ertheless, data do not represent a clear-cut textbook rendition of the RNGHI, 
but only an approximation. Therefore the bottom line of our combined theoret- 
ical and experimental research is this: The processes observed in EDFLs very 
likely constitute a manifestation of the RNGHI, but it is a manifestation in a 
'dressed' way. Inhomogeneous broadening, distributed losses, and the three- 
level structure of erbium ions take their imprint on the instability. Moreover, it 
is very likely that noise plays a central role in the 'intermittent' appearance of 
multimode emission. 

Theories are very often neater and much more elegant than real-world ex- 
periments. The closest thing to Lorenz-type laser chaos that was ever found 
experimentally suffered from complications that are absent from the model, as 
we commented in Sect. ^ In a similar way, the experiment in |89| presents the 
closest thing to RNGHI that has been found to date. In any case, we do not 
have any doubt that the experimental observations are a clear manifestation of 
the resonant Rabi instability, i.e., that Rabi sideband gain is the responsible for 
the observed instability. 

Another remarkable result of our research is that the RNGH mechanism 
is essential for understanding the multimode emission threshold in inhomoge- 
neously broadened lasers, even for relatively short values of the laser cavity 
depending on the active medium. This fact suggests that RNGHI could be im- 
portant for correctly understanding mode-locking in lasers with an accessible 
coherence length (see Eq. (|196|) and the subsequent discussion). 

Finally we would like to remark that there are open questions from both the 
theoretical and the experimental sides. On the one hand, the observed intermit- 
tent pulsations needs to be theoretically explained, and it must be determined 
up to what extent subcriticality and noise could explain them, or wether other 
phenomena we have not yet considered (such as dispersion or fiber nonlinear- 
ity |67[ I13p need to be taken into account. On the other hand, experimental 
research in other laser types, as NdDFLs or Nd:YAG lasers, would help to un- 
derstand how the RNGH mechanism affects real lasers. 
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